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FOREWORD 
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report  extends  the  nonlinear  AP95  code  to  the  roll  position  of  45  deg,  or  the  tail  fins  in  a  cross  (x) 
arrangement.  This  new  technology  will  be  transitioned  in  a  later  version  of  the  aeroprediction  code, 
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1.0  INTRODUCTION 


Many  of  the  world’s  missiles  fly  in  either  the  roll  stabilized  position  of  <E>  =  0  deg  (or  plus 
fin  orientation)  or  =  45  deg  (or  cross  fm  orientation).  Figure  1  illustrates  these  fin  orientations 
for  a  cruciform  missile  looking  from  the  missile  nose  toward  the  rear.  As  illustrated  in  the  figure, 
the  O  =  0  deg  plane  generally  gives  slightly  more  normal  force  and  a  slightly  more  stable 
configuration  in  pitch  at  a  given  angle  of  attack  than  does  the  missile  rolled  to  O  =  45  deg  (the 
physics  of  why  this  occurs  will  be  discussed  in  the  analysis  section).  On  the  other  hand,  a  missile 
in  the  <!>  =  45  deg  plane  is  in  a  roll-stable  position,  which  means  less  energy  is  required  to  maintain 
a  constant  roll  orientation.  Also,  all  four  fins  can  be  deflected  simultaneously,  giving  30  to 
50  percent  more  normal  force  from  control  deflection  than  only  two  fins  deflected  in  the  ^  =  0  deg 
roll  position. 

The  latest  version^  of  the  NSWCDD  Aeroprediction  Code  (ARC),  AP95,  calculates 
aerodynamics  only  in  the  <!>  =  0  deg  plane.  As  such,  aerodynamics  at  O  =  45  deg  roll  must  be 
obtained  from  another  aerodynamics  code^^  or  estimated  from  the  AP95  <D>  =  0  deg  results.  It  is  the 
intent  of  this  report  to  discuss  the  physical  phenomena  that  occur  at  the  ^  =  0  and  45  deg  roll 
orientations  and  to  develop  a  semiempirical  mathematical  model  that  will  allow  the  AP95  to  be 
extended  to  the  41>  =  45  deg  plane. 

The  AP95  is  an  approximate  analytical  Aeroprediction  Code  primarily  designed  to  provide 
preliminary  estimates  of  aerodynamics  for  use  in  particle  ballistic  models,  trim  aerodynamic  models, 
or  structures  and  heat  transfer  models.  Static  aerodynamics  are  generally  estimated  within 
±10  percent  and  center  of  pressure  within  ±4  percent  of  the  body  length.  To  obtain  aerodynamics 
estimation  accuracy  levels  generally  desired  for  full  six-degree-of-freedom  (6  DOF)  simulations 
requires  either  a  more  accurate  numerical  code^  or  wind  tunnel  data  or  both.  As  such,  the  AP95 
does  not  attempt  to  compute  out-of-the-pitch-plane  aerodynamics  or  coupling  effects  between  the 
pitch  and  yaw  planes.  On  the  other  hand,  if  one  is  interested  in  the  <1>  =  45  deg  plane,  and  control 
deflections  are  symmetric  with  respect  to  the  pitch  plane,  then  the  AP95  code  can  be  modified  to 
allow  aerodynamics  for  the  O  =  45  deg  plane.  These  modifications  will  be  made  in  an  analogous 
process  to  the  nonlinear  semiempirical  methodology  of  Reference  1 . 

One  may  ask,  why  not  use  another  engineering  code  such  as  that  of  Reference  2  or  3  to 
compute  aerodynamics  in  the  =  45  deg  plane.  There  are  several  reasons  for  this.  First  of  all,  it 
is  more  convenient  to  have  the  capability  for  O  =  0  and  45  deg  aerodynamics  in  a  single  code  as 
opposed  to  two  codes.  Second,  the  codes  of  References  2  and  3  are  more  limited  in  angle  of  attack 
and  Mach  number  than  the  AP95  code  at  41>  =  0  deg.  Third,  along  with  other  works,^  Reference  1 
has  shown  the  AP95  to  give  more  consistently  accurate  static  aerodynamics  in  the  =  0  deg  roll 
orientation  than  either  the  Reference  2  or  3  code.  It  is  believed  this  accuracy  level  can  be  extended 
to  the  ^>  =  45  deg  roll  orientation. 
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0=0° 


0=45° 


TYPICAL  FORCE  AND  MOMENT  COMPARISONS 


0=0° 


0=45° 


•  (Cn)o=o°  generally  higher  than 
(Cn)o=45°  at  higher  a 

•  GENERALLY  MORE  STABLE 
IN  PITCH  AT  0=0°  vs  0=45° 


•  NATURALLY  STABLE  POSITION 
IN  ROLL  (LESS  ENERGY  TO 
MAINTAIN  CONSTANT  ROLL) 

•  MORE  CONTROL  AUTHORITY 
AT  0=45°  DUE  TO  FOUR  FINS 
DEFLECTED  vs  TWO  AT  0=0° 


FIGURE  1.  CRUCIFORM  WING-BODY-TAIL  MISSILE  CONFIGURATION 
FLYING  AT  ROLL  OF  O  =  0  AND  45  DEG 
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On  the  other  hand,  the  methods  of  References  2  and  3  can  compute  aerodynamics  at  any  roll 
orientation.  Hence,  if  one  wants  to  perform  6  DOF  performance  analysis  and  is  willing  to  accept 
potentially  large  inaccuracy  levels  in  the  out-of-pitch-plane  aerodynamics  over  a  more  limited  range 
of  angle  of  attack  (AOA)  and  Mach  number,  then  References  2  and  3  are  the  only  engineering  codes 
available  in  the  United  States  for  performing  this  analysis,  to  the  authors’  knowledge.  The  approach 
of  the  AP95  is  to  neglect  the  body  vortex  scale  effects  due  to  length  and  leave  the  body  vortex 
effects  in  the  fm  aerodynamics.  This  eliminates  the  errors  associated  with  applying  a  semiempirical 
theory  to  subtract  out  nonlinear  body  vortex  effects  on  one  configuration  and  then  add  them  back 
in  for  a  different  configuration. 

The  body  vortex  effects  are  therefore  inherently  a  part  of  the  interference  factors  of  the  wing 
in  conjunction  with  the  body  and  the  body  in  conjunction  with  the  wing.  Based  on  limited 
comparisons  to  date,  the  direct  approach  of  AP95  of  using  data  from  wind  tunnel  results  versus 
modeling  data  for  the  wing  in  conjunction  with  the  body  for  body  vortex  effects  appears  to  give 
better  agreement  with  experiment  on  configurations  outside  the  wind  tunnel  data  base.  It  is  also 
much  more  straightforward  to  apply  and  requires  less  computational  time  as  well.  Certainly,  as  one 
attempts  to  achieve  angle-of-attack  capability  exceeding  20  deg,  the  errors  associated  with 
subtracting  out  body  vortex  effects  with  a  semiempirical  theory  based  on  limited  data  become  more 
and  more  questionable. 

The  overall  approach  to  modify  the  AP95  code  for  nonlinear  aerodynamics  in  the  =  45 
deg  plane  will  thus  be  very  similar  to  that  for  the  <E>  =  0  deg  plane.  Linearized  theories  or  slender 
body  theory  will  be  used  for  low  angle  of  attack  estimates  and  the  data  bases  of  References  6-9  will 
be  used  to  develop  empirical  or  semiempirical  corrections  to  account  for  the  nonlinearities  that  occur 
in  normal  force  and  center  of  pressure  with  increasing  angle  of  attack. 


2.0  ANALYSIS 


As  indicated  in  the  introduction,  the  goal  is  to  develop  a  nonlinear  semiempirical  model  for 
cruciform  missiles  for  the  <I>  =  45  deg  plane.  It  is  envisioned  that  this  model  will  be  analogous  to 
the  ^  =  0  deg  plane  methods  in  AP95,  except  the  normal  force  and  pitching  moments  due  to  the 
wing  alone  and  interference  aerodynamics  will  have  to  be  derived  for  the  3>  =  45  deg  roll 
orientation. 


Referring  to  the  total  normal  force  coefficient  equation  for  a  wing-body-tail  configuration 
as  given  by  Reference  9  we  have 


C  =  C  + 

N  '"Np 


(k  +K  +k 

W(B)  B(W)/  ^  W(B)  B(W)/  wj  \  N«/ 


W 


[(^T(B)  ^  ^B(T))  f'^T(B)  ^B(T))  ^x]  (‘^Noc)t 
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In  the  AP81/°  all  terms  in  Equation  (1)  were  linear  except  for  body-alone  normal  force.  As  a  result, 
the  AP81  gave  increasingly  erroneous  results  as  the  nonlinear  aerodynamics  associated  with  AOA 
became  important.  This  typically  meant  10  to  15  deg  AOA  as  an  upper  limit  The  AP93“’‘^ 
developed  nonlinear  corrections  for  each  of  the  terms  in  Equation  (1)  up  to  AOA  30  deg  using  wind 
tunnel  component  data  bases^'*.  A  more  consistently  accurate  Aeroprediction  Code  was  obtained 
up  to  AOA  30  deg.  Typical  average  accuracy  boundaries  of  ±10  percent  on  axial  and  normal  force 
and  ±4  percent  of  body  length  for  center  of  pressure  could  be  obtained.  The  AP95'  extended  the 
AOA  nonlinear  corrections  of  Equation  (1)  to  90  deg  AOA  through  extrapolation  of  the  data  bases 
of  References  6-8  plus  the  more  recent  data  base  of  Reference  9,  along  with  engineering  judgement. 
Jorgenson's*^  data  were  also  used  for  improvements  in  body-alone  aerodynamics.  The  accuracy 
levels  of  the  AP93  code  were  maintained  for  the  AP95  except  in  regions  of  strong  internal  shock 
interactions  (high  M,  high  AOA). 

The  way  Equation  (1)  is  implemented  in  both  the  AP93  and  AP95  codes  is  by  computing 
the  nonlinearities  of  the  body-alone  term  based  on  AOA  only,  but  all  the  other  nonlinear  corrections 
to  the  interference  factors  and  wing  or  tail  aerodynamics  are  based  on  the  total  angle  of  attack,  or 
«  ±  6.  This  was  done  for  ease  of  implementation  into  the  operational  AP81  code.  Hence,  even 
though  the  component  nonlinear  terms  were  derived  from  wind  tunnel  data  bases  from  AOA  data 
only,  they  were  implemented  in  a  total  local  AOA  sense.  The  empirical  nonlinear  corrections  were 
then  fine-tuned  based  on  comparisons  to  other  configuration  aerodynamics  outside  the  data  bases. 
The  center  of  pressure  of  individual  missile  component  aerodynamics  were  treated  similarly  to  the 
AP81  except  shifts  in  center  of  pressure  based  on  data  were  derived.  These  shifts  were  implemented 
in  tabular  form  as  a  function  primarily  of  AOA  and  M„. 

Also,  the  Cj^  term  of  Equation  (1),  which  is  the  downwash  normal  force  on  the  body  due 
to  the  wing  shed  vortices,  is  neglected.  This  is  because  it  is  inherently  included  in  the  wind  tunnel 
data  bases  and  it  is  believed  the  errors  in  trying  to  analytically  estimate  the  term,  subtract  it  out  on 
one  configuration,  and  then  add  it  back  in  later  on  a  different  configuration,  are  as  large  or  larger 
than  the  errors  from  incorporating  it  into  the  Kb(W)  term. 

Equation  (1)  can  also  be  rewritten  as 


Cn  = 


+  c 
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(2A) 


where  it  is  understood  that  encompasses  the  term.  For  ease  of  implementation  into  an 
existing  code  designed  primarily  for  linear  aerodynamics,  most  of  the  terms  in  Equation  (2A)  are 
separated  into  a  linear  and  nonlinear  contribution  due  to  «  or  6.  For  example,  the  wing-body  term 
is  computed  in  the  AP95  code  as  follows: 
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(2B) 


The  linear  or  small  angle  of  attack  terms  of  Equation  (2B)  are  estimated  by  linear  theory  (LT)  or 
slender  body  theory  (SBT).  This  gives  the  Aeroprediction  Code  a  good  fundamental  basis  for  its 
aerodynamic  estimates.  The  nonlinear  corrections  due  to  higher  AO  A  or  control  deflection  are  each 
estimated  directly  from  component  wind  tunnel  data  bases.®'^  Each  of  the  other  terms  in  Equation 
(2A)  is  treated  in  a  similar  fashion  to  Equation  (2B)  in  the  actual  implementation  into  the 
Aeroprediction  Code. 

In  the  context  of  Equation  (2),  we  therefore  seek  the  nonlinear  definition  of  each  of  the  terms 
in  Equation  (2).  It  is  expected  that  the  body-alone  term  [first  term  of  Equation  (2A)]  will  be 
independent  of  <I>.  In  reality,  this  is  not  necessarily  the  case  for  M<2  and  high  angle  of  attack 
because  of  the  asymmetric  shedding  of  vortices.  The  mechanism  of  this  shedding  is  not  clear,  but 
it  is  suspected  that  slight  imperfections  in  the  flow  or  body  shape,  from  uniform  or  axisymmetric, 
respectively,  could  contribute  to  this  phenomenon.  At  present,  the  Aeroprediction  Code  does  not 
account  for  out  of  plane  aerodynamics,  and  therefore  the  side  force  created  by  the  asymmetric 
shedding  of  body  vortices  is  not  predicted.  Also,  in  the  Reference  8  data,  normal  force  varied  by 
about  10  percent  as  a  function  of  roll  in  the  region  of  asymmetric  vortex  shedding.  Instead  of 
including  this  variation,  it  was  averaged  out. 

Each  of  the  remaining  terms  in  Equation  (2)  will  be  predicted  in  an  analogous  fashion  to  the 
AP95  code  developed  for  4>  =  0  deg,  except  here  the  quantities  will  be  for  O  =  45  deg.  As  already 
mentioned,  that  approach  was  based  on  linear  theory  or  slender  body  theory  for  small  values  of  « 
and  empirical  databases  (References  6-9)  to  develop  nonlinear  corrections  for  large  «.  As  such,  it 
is  instructive  to  examine  the  fundamental  impact  of  roll  orientation  on  linearized  and  slender  body 
theory  before  proceeding  to  the  nonlinear  corrections,  which  are  empirical  in  nature. 


2. 1  SLENDER  BODY  AND  LINEAR  THEORY  RESULTS  FOR 

ROLL-DEPENDENT  AERODYNAMICS 

References  14  and  15  were  primary  materials  used  for  examining  roll  dependence 
implications  from  slender  body  and  linearized  theories.  A  somewhat  detailed  summary  of  these 
results  is  given  for  information  purposes  in  Appendix  A.  The  summary  of  the  key  findings  in 
Appendix  A,  repeated  here  for  convenience,  are 

a)  For  cruciform  wings  alone  or  a  wing-body  combination,  the  total  normal  force  is 
independent  of  roll. 


5 


NSWCDD/TR-95/160 


b)  For  a  planar  wing-body  combination  at  roll,  the  loading  on  the  windward  plane  panel 
is  greater  by  an  equal  amount  to  that  on  the  leeward  plane  panel.  This  means  that 
if  one  were  trying  to  design  a  code  for  lateral  aerodynamics,  roll  dependence  of  each 
fm  planform  must  be  considered.  On  the  other  hand,  if  longitudinal  aerodynamics 
are  of  primary  interest,  the  total  normal  force  on  the  entire  wing  planform  can  be 
considered. 

c)  For  a  cruciform  wing-body-tail  configuration  at  roll,  eight  vortices  are  shed  in  the 
wing-body  region,  which  adversely  affects  the  tail  lift.  This  is  as  opposed  to  four 
vortices  at  <>  =  0  deg. 

d)  The  planar  theory  developed  for  wing-tail  interference  can  be  used  to  approximate 
the  loss  of  lift  on  the  tails  at  <J>  =  45  deg. 

e)  The  aerodynamics  of  a  cruciform  wing-body-tail  combination  with  zero  control 
deflections  are  independent  of  roll  position. 

These  findings  for  roll  dependence  from  linearized  or  slender  body  theory  are  quite  useful 
in  helping  plan  how  to  develop  a  nonlinear  APC  for  <I>  =  45  deg.  While  the  conclusions  of  linear 
theory  roll  dependence  may  not  translate  to  the  nonlinear  case,  we  will  still  use  the  findings  to  help 
guide  the  nonlinear  code  development.  In  particular,  the  item  (a)  conclusion  implies  use  of  the 
4>  =  0  deg,  wing  alone  data  for  <1^  =  45  deg.  This  is  quite  important  because  the  available  wing- 
alone  data  bases  are  all  at  =  0  deg.  This  means  that  any  nonlinear  wing-alone  roll  dependence 
will  be  included  in  the  interference  factors  rather  than  the  wing-alone  solution,  which  is  independent 


The  second  major  result  of  the  key  slender  body/linear  theory  roll  dependence  findings  is 
that  for  cruciform  missiles,  we  can  use  the  same  interference  approaches  as  in  the  AP95,  except  the 
constants  need  to  be  changed  because  of  a  different  roll  angle.  The  combination  of  these  two 
conclusions  are  quite  important  because  they  basically  allow  the  direct  usage  of  the  AP95  code  with 
different  constants  for  the  nonlinear  interference  terms  at  ^>  =  45  deg  versus  =  0  deg. 

The  third  significant  conclusion  is  that  for  small  AOA,  wing-body-tail  aerodynamics  are 
independent  of  roll  position.  This  allows  the  usage  of  wing-tail  interference  methodology  designed 
for  planar  computations  for  different  roll  orientations,  so  long  as  the  proper  number  of  vortices  are 
considered.  Again  different  nonlinear  corrections  are  expected  for  the  <I>  =  45  deg  versus  the 
<!>  =  0  deg  roll  position. 


2.2  NONLINEAR  AERODYNAMICS  METHODS 

This  section  will  describe  the  methods  used  for  computing  the  nonlinear  corrections  for  each 
of  the  terms  in  Equation  (2).  These  corrections,  with  the  exception  of  the  body  alone,  are  all 
empirical  in  nature. 
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2.2.1  Wing-Alone  Method 

One  of  the  primary  reasons  for  analyzing  the  slender  body  and  linear  theory  implications  on 
component  aerodynamics  in  Appendix  A  was  to  show  that  the  wing-alone  methodology  developed 
for  the  <I>  =  0  deg  plane  could  also  be  applied  in  the  «I>  =  45  deg  plane.  This  was  at  least  true  in  the 
linear  sense.  It  will  also  be  assumed  to  be  true  in  the  nonlinear  sense.  Any  nonlinearities  not 
accounted  for  in  using  the  <I>  =  0  deg  methods  for  <I>  =  45  deg  will  therefore  be  included  in  the 
interference  factors. 

As  a  result  of  this  approach,  the  methodology  of  Reference  16  can  be  used  directly.  This 
methodology  uses  a  fourth  order  equation  in  angle-of-attack  to  estimate  wing-alone  lift,  as  opposed 
to  the  second  order  approach  of  Reference  11.  The  specific  equations  are  defined  by 

+  “3%  +  ^4%  (3) 

>2  =  34.044  -  4.824  ^  0.426  -  6.412.,  (4) 

a,  =  -88.240  23.032  -  2.322(C„)__^^.  +  11.464.,  (5) 

..  =  53.219  -  17.595  .  2.661  -  5.971.,  (6) 


The  term  ai  of  Equations  (3)-(6)  is  the  value  of  the  wing-alone  lift  curve  slope  at «  =  0  given  by 
linear  theory.  The  terms  (Cjjj  (Cjj)  ,  -  „  are  values  of  the  wing-alone  normal  force 

coefficients  at  «=  15,  35,  and  60  deg,  respectively,  defined  by  the  data  bases  of  References  2,  6  and 
7.  Above  of  60  deg,  extrapolation  of  the  aerodynamics  at  of  60  deg  is  used.  For  more  details 
of  the  method,  the  reader  is  referred  to  Reference  16.  As  shown  in  Reference  16,  the  fourth-order 
method  of  Equations  (3)-(6)  improves  the  wing-alone  prediction  accuracy  below  =  30  deg  over 
the  second-order  method  used  in  the  AP93,  while  allowing  wing-alone  aerodynamic  estimation  to 
180  deg. 

The  value  of  Cj^  computed  by  Equation  (3)  includes  the  linear  and  nonlinear  term.  To 
include  it  into  a  term  like  Equation  (2B),  requires  this  term  to  be  separated  into  a  linear  and 
nonlinear  component.  This  is  easily  done,  as  the  linear  term  is  known;  so  the  nonlinear  term  is 
simply  the  difference  between  the  total  wing-alone  value  and  its  linear  counterpart.  The  secant  slope 
is  then  formed  for  each  of  these  terms  by  dividing  the  linear  and  nonlinear  parts  by  the  local  AO  A, 
«  +  6.  This  approach  was  taken  for  ease  of  incorporation  into  an  existing,  operational,  linearized 
code. 


The  center  of  pressure  of  wing-alone  normal  force  is  predicted  the  same  way  as  the  AP93. 
That  approach  used  a  second  order  method  to  move  the  center  of  pressure  from  the  linear  theory 
value  at «  =  0  to  the  centroid  of  the  planform  area  at «  =  60  deg. 
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22.2  Body- Alone  Method 

The  body-alone  methodology  is  the  same  as  that  of  Reference  1,  with  one  exception.  The 
one  exception  from  the  Reference  1  work  is  a  slight  shift  in  the  "drag  bucket"  at  low  crossflow 
Mach  numbers.  Reference  1  was  optimized  to  the  Reference  8  data  base,  which  seemed  to  show  a 
critical  crossflow  Reynolds  number  (the  point  where  the  C ^  curve  dropped  from  a  value  of  1 .2  to 
a  lower  value  that  was  caused  by  the  flow  remaining  laminar  over  a  greater  portion  of  the  body) 
around  R^  of  330,000.  On  the  other  hand,  most  other  wind  tunnel  results  seemed  to  imply  the 
Reference  13  drag  bucket  to  be  more  representative  of  most  data  than  that  assumed  as  the  baseline 
for  the  Aeroprediction  Code.  As  a  result,  this  drag  bucket  was  implemented  into  the  Aeroprediction 
Code.  The  effect  of  this  is  that  the  code  is  more  in  line  with  most  other  data  bases,  but  is  no  longer 
as  accurate  for  the  Reference  8  data  base.  The  user  still  has  the  option  of  shifting  the  drag  bucket, 
however.  Also,  it  was  found  that  when  two  sets  of  lifting  surfaces  were  present,  the  body  alone  data 
upon  which  this  "drag  bucket"  was  based  were  no  longer  valid,  and  a  of  1.2  at  low  subsonic 
crossflow  Mach  numbers  seemed  to  be  more  appropriate. 

Reference  1  basically  modified  the  Allen-Perkins*’  methodology  somewhat  by  using  sin  « 
versus  «  and  including  crossflow  Reynolds  number  in  the  crossflow  drag  computation  at  low 
crossflow  Mach  numbers.  The  center  of  pressure  estimation  is  the  same  as  that  used  in  the  AP95. 
In  that  approach,  the  nonlinear  normal  force  center  of  pressure  was  assumed  to  occur  at  the  centroid 
of  the  planform  area.  The  linear  theory  center  of  pressure  was  computed  analytically  in  most  cases. 
Finally,  an  empirical  data  base  was  used  to  approximate  center  of  pressure  shifts  due  to  transonic 
and  asymmetric  vortex  effects. 


2.2.3  Wing-Bodv  and  Bodv-Wine  Interference  Due  to  AO  A 


The  wing-body  and  body-wing  interference  factors  were  computed  using  a  combination  of 
the  Reference  8  data  base  for  0.6  ^  M„  ^  4.6  and  the  Reference  9  data  base  at  M  =  0.1.  Outside 
these  Mach  limits,  extrapolations  were  made  to  allow  the  methodology  to  compute  aerodynamics 
at  all  Mach  numbers.  These  extrapolations  were  not  as  difficult  as  they  may  seem  because  normal 
forces  and  center  of  pressure  have  basically  leveled  out  at  M  =  4.6  and  further  increases  in  Mach 
number  produce  fairly  small  changes  in  these  parameters. 


The  wing-body  interference  factor  is  defined  as 

C. 
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Here  Cj^  was  measured  directly  in  the  data  base  of  Reference  8  by  having  the  wing  in  close 
proximity  to  the  body  and  measuring  directly  the  load  on  the  wing  in  the  presence  of  the  body. 
Since  the  normal  force  was  measured  normal  to  a  single  fm,  to  get  the  normal  force  on  the  wing  in 
the  presence  of  the  body  at  =  45  deg  from  the  data  required  the  data  to  be  multiplied  by  cos  «I>. 
To  reduce  measurement  errors,  the  data  from  all  four  fins  were  averaged.  Hence,  if  subscripts  1 
through  4  represent  the  loads  on  the  individual  fins,  then 


cos 

2 


(8) 


No  attempt  was  made  to  correct  for  wind  tunnel  errors  near  zero  AOA  caused  by  flow 
misalignments.  These  errors  can  cause  the  normal  force  curve  to  be  shifted  as  much  as  a  degree. 
This  means  that  the  accuracy  could  have  some  slight  errors  near  zero  AOA. 


Based  on  the  accuracy  analysis  of  Reference  1,  this  should  give  fairly  accurate  values  of 
Kw(B)  for  all  but  the  highest  aspect  ratio  where  the  wing  planform  area  was  only  about  2  percent  of 
the  body  planform  area  in  the  crossflow  plane.  Cjj  of  Equation  (7)  was  arrived  at  from 
Reference  1,  which  in  turn  used  the  data  bases  of  References  2,  6  and  7.  Figure  2  shows  the  wing 
and  wing-body  configurations  tested  in  References  6,  8,  and  9,  respectively. 
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FIGURE  2A.  MODELS  USED  IN  LANGLEY®  WING-BODY  TESTS 
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FIGURE  2B.  MODELS  USED  IN  STALLINGS  AND  LAMB®  WING-ALONE  TESTS 


The  body-wing  interference  factor  is  defined  as 


K 


B(W) 


"N, 


B(W> 


(9) 


Unfortunately,  was  not  a  quantity  that  was  measured  directly  in  Reference  6  but  was 
computed  from  three  other  independent  measurements  of  body  alone,  wing  in  conjunction  with  the 
body  and  total  normal  force.  The  computation  for  was  then  made  by 

(10) 

As  shown  in  the  Reference  1  error  analysis,  this  process  gave  potential  errors  that  were  much  higher 
than  for  C„  ,  particularly  for  the  smaller  wings  (AR  ^  1 .0)  and  higher  Mach  numbers  (M  ^  2.5), 

^W(B) 

where  the  term  decreased  to  the  point  it  was  within  the  accuracy  of  the  data.  As  a  result, 
much  more  scatter  in  the  data  is  expected  for  this  term  and  more  engineering  judgement  is  required 
in  the  empirical  model  development. 

Unlike  Reference  8,  Reference  9  had  a  fairly  large  wing  planform  area  compared  to  the  body 
planform  area  (approximately  60  percent).  Moreover,  was  apparently  measured  separately. 

^B(W) 

Hence,  C.,  could  be  computed  based  on  direct  measurements  as  well  as  calculated  similar  to 
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Equation  (10).  Also,  data  were  obtained  all  the  way  to  «  =  90  deg  and  at  an  r/s  value  of  0.25.  As 
a  result,  more  confidence  is  placed  on  the  Reference  9  body-wing  interference  at  high  angle  of 
attack  than  the  Reference  8  data.  Unfortunately,  the  Reference  9  data  were  taken  only  at  M„  =  0.1, 
so  it  is  hard  to  extrapolate  it  past  about  =  0.6.  Fortunately,  it  compliments  the  Mach  number 
range  of  the  larger  Reference  8  data  base  quite  nicely. 

Figures  3-21  present  the  computed  results  for  Kw(b)  and  Kb(w)  for  the  various  fm  planforms 
at  the  various  Mach  numbers  covered  by  the  data  bases  of  References  8  and  9  for  ^  =  45  deg.  Also 
shown  on  the  figures  are  the  slender  body  theory  or  linear  theory  results  for  Kw(b)  and  Kb(W)  along 
with  approximations  to  the  data  that  will  be  used  in  the  empirical  model.  In  examining  the  ^>  =  45 
deg  data  of  Figures  3-21,  it  is  seen  that  the  data  follow  the  same  general  trends  as  the  «I>  =  0  deg  data 
of  Reference  1.  That  is,  Kw(b  decreases  with  increasing  angle  of  attack  until  a  minimum  is  reached. 
The  Reference  9  data  showed  this  minimum  to  be  about  20  percent  lower  than  the  «I>  =  0  deg  data. 
This  difference  in  minimums  could  not  be  confirmed  by  the  Reference  6  data  because  of  a  lack  of 
data  above  ec  =  25  to  30  deg  subsonically  and  40  deg  supersonically.  In  fact,  at  supersonic  speeds, 
the  Reference  6  data  seemed  to  indicate  a  minimum  value  of  Kw(B)  near  one  at  high  «,  similar  to  the 
<E>  =  0  deg  data.  Also  as  seen  in  the  figures,  SBT  gives  a  reasonable  approximation  to  Kw(b)  at «  up 
to  about  10  to  15  deg,  but  becomes  increasingly  inaccurate  as  «  increases. 

The  Kb(w)  data  follow  generally  the  same  trends  as  that  for  =  0.  However,  it  is  seen  that 
the  magnitude  of  Kb(w)  is  somewhat  different  at  both  subsonic  and  supersonic  speeds.  Take  a  couple 
of  specific  cases  for  example  to  illustrate  this.  One  of  these  is  the  M  =  0.8  data  for  AR  =  0.5.  The 
=  0  deg  data  increase  with  «  until  a  maximum  is  reached  somewhere  around  20  to  30  deg  AO  A, 
after  which  Kb(W)  starts  decreasing.  The  maximum  value  is  above  that  estimated  by  SBT.  *  Referring 
to  Figure  5  for  the  ^>  =  45  deg  data,  it  is  seen  that  the  Kb(w)  data  start  out  near  the  SBT  but  decrease 
continually  with  angle  of  attack.  The  minimum  value  of  25  percent  of  SBT  is  assumed  based  on 
the  Figure  3  high  AOA  data.  This  fairly  large  difference  between  the  =  0  and  ^>  =  45  deg  data 
for  Kb(W)  is  the  main  reason  missiles  flying  at  the  «1>  =  45  deg  orientation  have  less  lift  and  static 
stability  than  those  flying  at  ^>  =  0  deg  orientation  at  subsonic  Mach  numbers  and  ■>'  >  10  deg. 

A  second  case  is  considered  to  compare  Kb(W)  at  =  0  and  =  45  deg.  Again,  the  aspect 
ratio  0.5  case  is  considered,  but  at  M„  =  3.0  (see  Figure  10).  Reference  1  gave  the  3>  =  0  deg  Kb(W) 
as  close  to  the  SBT  value  at  «  =  0,  decreasing  to  0.4  to  0.6  at  «  =  25  deg.  Again  Figure  10  shows 
a  value  of  Kb(W)  of  about  half  that  at  ^>  =  0  deg,  again  indicating  a  lower  normal  force  and  less 
stability  from  the  body  carryover  term  for  $  45  deg  versus  $  =  0  deg  roll. 

All  the  data  in  Figures  3-21  for  Kw(B)  and  Kb(W)  can  be  generically  represented  by  Figures 
22  and  23.  Referring  to  Figure  22  first,  it  is  seen  that  there  are  six  parameters  needed  to  represent 
the  Kw(B)  as  a  function  of  «,  M,  and  fm  planform.  These  are  defined  as  follows: 

^  =  difference  between  SBT  and  data  at  <='  =  0 

=  angle  of  attack  where  Kw(b)  starts  decreasing 
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FIGURE  3.  WING-BODY  AND  BODY- WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  0. 1,  r/s  =  0.25,  AR  =  1 .5) 
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FIGURE  4.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  0.6,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  5.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  0.8,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  6.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (TVL  =  1-2,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  7.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  1.5,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  8.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  2.0,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  9.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (IVL  =  2.5,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  10.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  3.0,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  1 1 .  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  3 .5,  r/s  =  0.5,  AR  =  0.5) 
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FIGURE  13 .  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M.  =  0.6,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  14.  WING-BODY  AND  BODY-WING  EMTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  0.8,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  15.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (IvL  =  1.2,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  16.  WING-BODY  AND  BODY- WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  1.5,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  17.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M.  =  2.0,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  18.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  2.5,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  19.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M«,  =  3.0,  r/s  =  0.5,  AR=  1,2) 
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FIGURE  20.  WING-BODY  AND  BODY- WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M„  =  3.5,  r/s  =  0.5,  AR  =  1,2) 


NSWCDD/TR-95/160 


31 


FIGURE  21.  WING-BODY  AND  BODY-WING  INTERFERENCE  FACTORS 
AS  A  FUNCTION  OF  AOA  (M.  =  4.5,  r/s  =  0.5,  AR  =  1,2) 
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FIGURE  22.  GENERIC  REPRESENTATION  OF  WING-BODY  INTERFERENCE  FACTOR  AT  $  =  45  DEG 
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FIGURE  23 .  GENERIC  REPRESENTATION  OF  BODY- WING  INTERFERENCE  FACTOR  AT  $  =  45  DEG 
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=  rate  of  decrease  of  Kw(b)  between  «  =  «c  «  =  «d 
=  angle  of  attack  where  Kw(B)  reaches  an  initial  minimum 

=  angle  of  attack  where  Kw(b)  reaches  a  constant  value 


[^w(B)L  value  of  K^cb)  where  a  constant  value  with  angle  of  attack  has  been  attained 

M 

Tables  1-6  give  values  of  these  parameters  based  on  the  approximations  in  Figures  3-21. 


Figure  23  shows  that  Kb(W)  can  be  approximated  by  five  parameters.  These  are  defined  as 
follows: 


^^B(W) 

d- 


=  difference  between  SBT/LT  and  data  at «  =  0 
=  rate  of  change  of  Kb(W)  between  »<:  =  0  and  «  = 

=  angle  of  attack  where  dKg^j/d«  changes  sign 
=  angle  of  attack  where  Kb(W)  reaches  a  constant 

=  constant  value  of  Kb(w)  above  «  =  «2  as  a  percent  of  linear  theory  or  slender  body 
theory  (see  Figure  24) 


Tables  7-11  give  values  of  these  parameters  as  a  function  of  M„  and  the  fin  planform  parameters. 


The  mathematical  models  for  Kw(b)  and  Kb(W)  are  once  again  defined  based  on  SBT/LT  and 
using  the  data  in  Tables  1-11.  The  specific  equations  for  Kw(b)  are 
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TABLE  1,  DATA  FOR  ^  AT  <I»  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^5.0 

^0.25 

0, 0.5, 1.0 

-0.20 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

0.00 

0.05 

0.00 

-0.13 

0.00 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

0.13 

0.05 

0.00 

0.00 

-0.10 

0 

0 

0 

0 

0 

0 

^2.0 

0.5 

0.00 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

0.5 

0 

0.00 

0.10 

0.14 

0.00 

-0.22 

0 

0 

0 

0 

0 

0 

^2.0 

0 

0.00 

0.00 

0.00 

0.00 

-0.18 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

0.00 

0.10 

0.05 

-0.23 

0.00 

0 

0 

0 

0 

0 

0 

^2.0 

1.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0 

0 

0 

0 

0 

0 

TABLE  2.  DATA  FOR  AT  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^5.0 

^0.25 

0,0.5,  1.0 

10.0 

22.0 

22.0 

0.0 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

45.0 

11.5 

11.0 

10.0 

0.0 

0 

0 

0 

0 

0 

0 

1.0 

0.5 

45.0 

13.3 

0.0 

6.5 

0.0 

0 

0 

0 

0 

0 

0 

k2.0 

0.5 

20.0 

10.0 

0.0 

6.5 

2.2 

0 

0 

0 

0 

0 

0 

0.5 

0 

39.0 

15.0 

11.5 

10.0 

0.0 

0 

0 

0 

0 

0 

0 

^2.0 

0 

20.0 

10.0 

0.0 

6.5 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

45.0 

15.0 

15.0 

15.0 

0.0 

0 

0 

0 

0 

0 

0 

^2.0 

1.0 

20.0 

10.0 

0.0 

6.5 

1.5 

0 

0 

0 

0 

0 

0 
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TABLE  3 .  DATA  FOR  AT  $  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^0.25 

0, 0.5, 1.0 

1.0 

1.0 

1.00 

LOO 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

0.5 

1.0 

1.0 

1.00 

0.90 

0.90 

1.00 

0.95 

1.00 

0.97 

1.00 

1.0 

1.0 

0.5 

1.0 

1.0 

1.00 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

^2.0 

0.5 

1.0 

1.0 

0.95 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

0 

1.0 

1.0 

1.00 

1.00 

0.90 

0.90 

0.90 

0.90 

0.90 

0.90 

1.0 

^2.0 

0 

1.0 

1.0 

0.95 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

0.5 

1.0 

1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

^2.0 

1.0 

1.0 

1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

0.93 

0.90 

0.95 

1.0 

TABLE  4.  DATA  FOR  AT  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

sO.l 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^0.25 

0,0  .5, 1.0 

20.0 

40.0 

38.0 

35.0 

30.0 

25.0 

16.3 

15.1 

13.9 

13.1 

10.0 

0.5 

0.5 

59.0 

33.0 

30.0 

25.6 

25.0 

15.0 

15.0 

10.0 

15.0 

12.0 

10.0 

1.0 

0.5 

59.0 

38.0 

32.0 

26.0 

24.0 

17.0 

15.0 

14.4 

10.0 

10.0 

10.0 

^2.0 

0.5 

39.0 

31.5 

30.0 

28.0 

25.0 

16.5 

15.0 

14.4 

10.0 

13.0 

10.0 

0.5 

0 

39.0 

35.5 

33.0 

39.5 

29.5 

15.0 

25.0 

15.0 

15.0 

10.0 

10.0 

s:2.0 

0 

39.0 

31.5 

30.0 

28.0 

24.7 

17.0 

13.5 

11.4 

10.0 

10.0 

10.0 

0.5 

1.0 

59.0 

35.5 

33.0 

25.6 

29.5 

15.0 

15.0 

15.0 

12.0 

13.0 

10.0 

^2.0 

1.0 

39.0 

31.5 

30.0 

28.0 

23.3 

14.0 

16.0 

15.0 

11.8 

12.0 

10.0 
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TABLE  5.  DATA  FOR  AT  4*  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^0.25 

0,0.5, 1.0 

40.0 

45.0 

45.0 

40.0 

44.0 

43.0 

38.0 

28.0 

25.0 

29.0 

20.0 

0.5 

0.5 

65.0 

33.0 

30.0 

49.0 

52.0 

40.0 

40.0 

30.0 

25.0 

25.0 

20.0 

1.0 

0.5 

65.0 

38.0 

47.0 

49.5 

66.0 

48.5 

45.0 

41.0 

40.0 

10.0 

20.0 

^2.0 

0.5 

40.0 

31.5 

40.0 

56.0 

57.0 

45.0 

45.0 

41.0 

40.0 

28.0 

20.0 

0.5 

0 

40.0 

35.5 

33.0 

65.0 

48.0 

50.0 

46.0 

30.0 

30.0 

50.0 

20.0 

^2.0 

0 

40.0 

31.5 

40.0 

56.0 

55.0 

58.5 

49.8 

44.2 

41.5 

28.5 

20.0 

0.5 

1.0 

65.0 

35.5 

33.0 

49.0 

52.0 

40.0 

28.0 

24.0 

21.0 

13.0 

20.0 

^2.0 

1.0 

40.0 

31.5 

40.0 

56.0 

49.5 

44.0 

40.0 

33.0 

32.0 

28.0 

20.0 

TABLE  6.  DATA  FOR 


AT4>  =  45  DEG 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^0.25 

0,0.5,  1.0 

0.80 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

.5 

0.80 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.5 

0.80 

0.90 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

^2.0 

.5 

0.80 

0.90 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

0 

0.85 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

^2.0 

0 

0.85 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.5 

1.0 

0.80 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

^2.0 

1.0 

0.80 

0.95 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 
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TABLE?.  DATA  FOR  ^ 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

B 

^6.0 

0, 0.5, 1.0 

-0.10 

-0.18 

0.00 

0.00 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

0.5 

0.0 

-0.22 

0.00 

0.00 

0.2 

0.1 

0.08 

0 

0 

0 

0 

1.0 

0.5 

0.0 

-0.07 

-0.18 

0.20 

0.2 

0 

0 

0 

0 

0 

0 

^2.0 

0.5 

0.0 

-0.23 

-0.18 

0.20 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

0 

0.0 

-0.12 

0.00 

0.25 

0.2 

0 

0 

0 

0 

0 

0 

^2.0 

0 

0.0 

-0.23 

-0.18 

0.20 

0.0 

0 

0 

0 

0 

0 

0 

0.5 

1.0 

0.0 

-0.22 

0.00 

0.25 

0.2 

0.1 

0.08 

0 

0 

0 

0 

^2.0 

1.0 

0.0 

-0.23 

-0.18 

0.20 

0.0 

0 

0 

0 

0 

0 

0 

Here, 


d- 


where 

The  equations  for  the  body-wing  interference  factor  are 


K, 


BCW) 


[^bwLin}  ’  "l  ^  -  "2 


K  =  [K  1  •  OC  >  <x 

“‘BCW)  L  >  -  2 


(12) 
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TABLE  9.  DATA  FOR  (DEG) 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^.025 

0,0.5, 1.0 

10.0 

45.5 

35.0 

30.0 

23.0 

22.0 

20.8 

20.0 

20.0 

15.0 

8.0 

0.5 

0.5 

10.0 

57.0 

45.0 

30.0 

25.0 

16.0 

20.0 

15.0 

10.0 

10.0 

8.0 

1.0 

0.5 

10.0 

34.5 

35.0 

35.0 

36.6 

10.0 

15.0 

17.5 

42.0 

40.0 

30.0 

^2.0 

0.5 

10.0 

15.0 

30.0 

35.0 

20.0 

20.0 

18.0 

17.5 

30.0 

35.0 

35.0 

0.5 

0 

10.0 

35.0 

45.0 

30.0 

19.0 

20.0 

22.5 

15.0 

15.0 

15.0 

8.0 

^2.0 

0 

20.0 

15.0 

30.0 

35.0 

20.0 

20.0 

18.0 

17.5 

30.0 

35.0 

35.0 

0.5 

1.0 

10.0 

57.0 

45.0 

30.0 

25.0 

16.0 

20.0 

15.0 

10.0 

10.0 

8.0 

^2.0 

1.0 

10.0 

15.0 

30.0 

35.0 

20.0 

10.0 

15.0 

17.5 

30.0 

35.0 

35.0 

TABLE  10.  DATA  FOR  (DEG) 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

■ 

^6.0 

^0.25 

0, 0.5,  1.0 

55.0 

55.0 

50.0 

50.0 

45.0 

37.0 

44.0 

29.5 

29.5 

29.5 

25.0 

0.5 

0.5 

75.0 

65.0 

55.0 

43.0 

40.0 

38.0 

44.0 

44.0 

36.0 

30.0 

20.0 

1.0 

0.5 

75.0 

65.0 

60.0 

60.0 

60.0 

60.0 

62.0 

80.0 

42.0 

40.0 

30.0 

2:2.0 

0.5 

75.0 

65.0 

60.0 

60.0 

50.0 

60.0 

62.0 

80.0 

42.0 

45.0 

45.0 

0.5 

0 

75.0 

60.0 

60.0 

52.0 

40.0 

35.0 

44.0 

50.0 

36.0 

30.0 

20.0 

^2.0 

0 

75.0 

65.0 

60.0 

60.0 

50.0 

60.0 

62.0 

80.0 

42.0 

45.0 

45.0 

0.5 

1.0 

75.0 

65.0 

55.0 

42.0 

40.0 

38.0 

44.0 

40.0 

36.0 

30.0 

20.0 

^2.0 

1.0 

75.0 

65.0 

60.0 

60.0 

60.0 

60.0 

62.0 

80.0 

42.0 

45.0 

45.0 
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TABLE  11.  DATAF0R[K3^j|^  (FRACTION  OF  SBT/LT) 


MACH  NUMBER 

ASPECT 

RATIO 

TAPER 

RATIO 

^0.1 

0.6 

0.8 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

B 

^6.0 

^0.25 

0, 0.5, 1.0 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

0.5 

0.5 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

1.0 

0.5 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

^2.0 

0.5 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

0.5 

0 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

^2.0 

0 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

0.5 

1.0 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

S;2.0 

1.0 

0.25 

0.25 

0.25 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0 

Linear  interpolation  is  used  in  Tables  1-1 1  to  define  the  parameters  in  Equations  (11)  and 
(12)  other  than  [k^(b)1  and  ,  as  a  function  of  ML,  AR,  and  A..  These  two  terms  are 

calculated  analytically  using  SBT  or  LT  as  done  before  in  the  previous  versions  of  the 
Aeroprediction  Code. 

The  methodology  of  the  O  =  45  deg  technology  is  similar  to  that  for  <&  =  0  deg  (see 
Reference  1);  the  difference  is  in  all  the  constants  of  Tables  1-11.  Exceptions  to  this  are  "fm 
choking,"  r/s  effects,  and  center  of  pressure  shifts.  These  will  be  discussed  in  the  following  three 
sections. 


2.2.4  Nonlinearities  Due  to  Internal  Shocks 


Three  nonlinearities  are  difficult  to  account  for  with  available  data.  These  phenomena  occur 
at  moderate  (M>2.0)  and  higher  supersonic  Mach  numbers  and  higher  values  of  AOA  (typically 
^  25  deg).  The  first  nonlinearity  is  sometimes  referred  to  as  "fin  choking";  the  second  is  when  the 
bow  shock  intersects  the  fins;  and  the  third  is  shocks  from  forward-located  fins  intersecting  aft- 
mounted  fins.  Each  of  these  phenomena  will  be  discussed  individually.  "Fin  choking"  is  a 
phenomenon  similar  to  that  which  occurs  when  an  inlet  becomes  unstarted  or  a  wind  tunnel  achieves 
its  maximum  rate  of  flow  (an  increase  in  power  produces  no  more  mass  flow  through  the  inlet).  As 
the  body  increases  in  angle  of  attack  with  the  fins  oriented  in  the  X  or  cross  orientation,  the  flow 
between  the  fins  will  eventually  "choke"  at  some  angle  of  attack  and  at  moderate  to  large  supersonic 
Mach  numbers.  When  this  happens,  a  strong  shock  is  formed  just  in  front  of  the  producing 

a  high  pressure  region  on  the  fins  and  body.  This  high  pressure  region  is  shifted  forward  from 
where  it  would  be  if  supersonic  flow  occurred  through  the  fins.  While  the  absolute  value  of  pressure 
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on  the  body  is  higher  than  for  the  unchoked  flow,  it  occurs  over  a  much  smaller  region  and  hence 
gives  only  slightly  higher  body-wing  interference  lift. 

Modified  Newtonian  Theory  (MNT)  can  be  used  to  help  approximate  the  combination  of  M 
and  «  that  will  cause  fin  choking  to  occur  in  the  $  45  deg  plane.  MNT  is  defined  by“ 

Cp  =  sin^  (6^^  (13) 


where 


2  . 

(Y+1)  M." 

Y 

Y-l 

Y+1 

1 

Y-l  _1 

YMi 

2 

2ymJ  -  (Y-l) 

(14A) 


6^^  =  sin’^  (sin  0  cos  «  -  sin  «  cos  O  cos  0) 


(14B) 


and  0  is  the  local  slope  of  the  body  or  fin  surface  with  respect  to  the  body  axis.  The  local  surface 
Mach  number  can  be  approximated  by 


Ml  =  < 


/  2  ^ 


Y-i; 


Y-l 

Y 


-1 


L/ 


1/2 


where 


(15) 


+  1 


(16) 


Terms  Cp  and  Cp^  of  Equation  (16)  come  from  Equations  (13)  and  (14A),  respectively. 

Figure  25  gives  the  local  Mach  number  results  just  upstream  of  the  fins  using  Equation  (15). 
Three  Mach  numbers  were  considered  in  this  figure:  the  Mach  number  of  the  Reference  19  data  of 
2.7,  plus  a  lower  value  of  1.5,  and  an  upper  value  of  ». 
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It  is  interesting  to  note  that  as  AOA  increases,  freestream  Mach  number  has  less  and  less  of 
an  effect  on  local  Mach  number.  Also  note  that  the  flow  in  the  windward  plane  goes  subsonic,  even 
with  no  fins  present,  at  =  35,  43  and  47  deg  for  the  Mach  numbers  M„  =  1.5,  2.7  and  <», 
respectively.  This  low  Mach  number  flow  will  then  intersect  the  fins.  If  the  fins  are  in  the  <>=  1 3  5 
and  225  deg  planes,  internal  shocks  can  reflect  off  opposite  fins  if  r/s  is  small.  This  lowers  the  Mach 
number  further.  The  larger  the  fins,  the  more  reflections  and  more  blockage  that  will  occur  and 
hence  the  lower  the  value  of  M„  and  «  required  for  "fin  choking."  Also,  the  larger  the  fins,  the 
harder  it  is  for  the  flow  in  the  region  of  the  fins  to  escape  around  the  fins. 

Reference  18  actually  showed  a  stagnation  pressure  occurring  on  the  windward  side  of  the 
fins  at  =  135  and  225  deg  at «  =  35  deg  and  M„  =  2.7.  It  is  suspected  that  the  flow  on  the  body 
attempts  to  wrap  around  the  body  in  the  vicinity  of  the  fins  but  is  constrained  by  the  fins  and  the 
flow  coming  in  the  opposite  direction  down  the  fins.  These  two  flows  then  merge  and  cause  a  local 
stagnation  region  (which  means  zero  velocity).  The  schlieren  photograph  for  this  condition  showed 
a  strong  shock  ahead  of  the  fins  merging  with  the  bow  shock.  When  the  shock  ahead  of  the  fins 
merged  with  the  bow  shock,  the  bow  shock  was  pushed  further  away  from  the  body. 

It  is  suspected  that  fin  size  plays  a  significant  role,  along  with  angle  of  attack  and  Mach 
number,  in  "fin  choking."  This  statement  is  based  on  the  fact  that  wings  with  smaller  r/s  values 
(wing-dominated  configurations)  seem  to  have  more  nonlinear  aerodynamics  in  the  "fin  choking" 
range  of  M  and  «  than  do  those  configurations  that  have  larger  values  of  r/s  (body  dominated). 

Two  approximations  are  being  made  to  take  into  account  some  of  the  nonlinearities 
associated  with  high  angle-of-attack  aerodynamics.  First,  the  criteria  for  determining  when  to  use 
slender  body  or  linear  theory  are  strictly  used  only  for  AOA  less  than  15  deg.  That  is,  if 


P  AR  (1  +  A)  J-  +  1  >  4 

mp  ; 

and  K  ^  15  deg,  then  linear  theory  is  used  for  defining  [Kb(w)]min-  If  the  above  conditions  are  not 
satisfied,  [Kb(W)]min  is  determined  as  a  fraction  of  slender  body  theory.  This  has  the  overall  effect 
of  slightly  increasing  the  minimum  value  of  body-wing  lift  for  lower  values  of  r/s. 

The  second  change  is  to  estimate  the  region  of  "fin  choking"  based  on  MNT  and  then  use 
this  region  along  with  Reference  8  and  18  results  to  model  the  r/s  effect.  Figure  26  gives  the 
approximate  boundary  for  fin  choking  using  MNT  and  the  region  over  which  it  appears  the 
nonlinear  effect  for  Kb(w)  for  r/s  ^  0.2  is  defined.  The  specific  approach  for  implementation  into 
the  APC  is  discussed  in  the  section  on  r/s  effects  (2.2.5). 

The  second  nonlinear  problem  mentioned  is  the  penetration  of  the  bow  shock  by  the 
windward  plane  fins.  This  occurs  for  the  fin  in  the  =  1 80  deg  plane  when  in  the  plus  or  <!>  =  0  deg 
roll  orientation  before  it  does  when  the  vehicle  is  rolled  to  a  «I>  =  45  deg  position.  The  impact  of  this 
phenomenon  on  the  aerodynamics  is  not  clear  in  examining  Tables  1-11.  In  examining  fin  pressure 
data  for  the  two  cases  of  Reference  18,  it  was  found  that  the  portion  of  the  fin  within  the  shock  layer 
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FIGURE  26.  Kw(b)  AND  Kb(w)  PREDICTION  FOR  ^45  DEG  AND  r/s  ^0.2 
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had  a  much  higher  pressure  than  the  part  outside  the  shock  layer.  This  certainly  would  imply  a 
difference  in  aerodynamics,  and  in  fact  this  difference  is  certainly  a  part  of  Tables  1-11.  However, 
it  is  not  clear  from  the  tables  the  contribution  of  this  phenomenon. 

The  third  internal  shock  problem  is  the  impact  that  a  shock  from  a  forward  located  fm  has 
on  a  rear  mounted  fm.  This  phenomenon,  unlike  the  previous  two,  is  not  included  in  the  data  in 
Tables  1-11.  This  is  because  the  data  sets  in  both  References  8  and  18  had  only  one  set  of  fins 
present  and  hence  this  phenomenon  cannot  be  accounted  for. 

In  summary,  it  is  fair  to  say  that  while  the  present  approximate  method  takes  into  account 
some  of  the  nonlinearities  in  missile  aerodynamics  due  to  internal  shocks,  it  is  far  from  complete. 
Additional  new  data,  combined  with  numerical  computations,  are  needed  to  develop  an  approximate 
engineering  model  to  represent  these  nonlinearities  accurately  and  hence  completely  cover  the 
desired  Mach  number  and  angle  of  attack  boundary  of  interest  for  all  configurations  of  interest  to 
weapons  designers. 


2.2.5  Treatment  of  Nonlinearities  Due  to  r/s 

Figures  3-21  and  Tables  1-1 1  were  based  mostly  on  the  data  of  Reference  8  that  were  for  r/s 
of  0.5.  This  data  base  showed  little  identifiable  effect  from  the  "fm  choking"  and  other  phenomena 
shown  by  the  Reference  18  data  base,  which  was  based  on  a  r/s  of  0.2.  Based  on  this,  an  hypothesis 
is  assumed  that  the  "fm  choking"  phenomenon  is  more  important  for  wing  dominated  configurations 
where  r/s  is  small.  To  handle  these  situations,  the  following  strategy  was  followed  in  treating  the 
effects  of  r/s.  Hopefully,  this  strategy  will  be  validated  or  improved  upon  with  additional  wind 
tunnel  or  computational  results  at  values  of  r/s  less  than  0.5. 

r/s  >  0.5.  For  this  range  of  r/s.  Figures  3-21  and  Tables  1-1 1  are  used  directly  to  develop 
the  nonlinear  corrections  to  slender  body  or  linear  theory.  The  data  base*  is  based  on  an  r/s  of  0.5 
but  will  be  assumed  to  apply  for  r/s  >  0.5  as  well,  unless  validation  against  cases  outside  the  range 
of  the  data  base  proves  unacceptable. 

r/s  ^  0.2.  When  r/s  =  0.2,  the  limited  data  base  of  Stallings,  et  al.^*  is  used  along  with  MNT 
to  determine  if  "fm  choking"  or  fm-fm  interference  effects  are  present  or  not.  This  is  based  on 
Figure  26.  If  these  effects  are  present,  then  the  Kw(b)  and  Kb(W)  for  =  0,  r/s  ^  0.2  condition  is  used 
for  <I>  =  45  deg  up  to  the  AOA  where  fm  choking  is  important.  This  is  based  on  the  fact  that  for  low 
r/s  values,  the  total  wing-body  normal  force  appeared  to  be  basically  the  same  for  3>  =  0  and 
<I>  =  45  deg  for  the  AOA  below  where  "fm  choking"  or  fm-fm  interference  occurs.  While  the  data 
used  are  for  r/s  =  0.2,  the  assumption  will  be  made  that  these  data  hold  for  r/s  <  0.2  as  well. 

If  "fm  choking"  does  occur,  the  Kw(b)  and  Kb(W)  for  «I>  =  0  is  used  up  to  the  lower  boundary 
on  Figure  26  where  "fm  choking"  becomes  important.  For  AOAs  above  the  upper  boundary  on 
Figure  26,  the  <I>  =  45  deg,  r/s  ^0.5  methodology  is  used  for  Kw(B)- 
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The  Kb(w),  ^>  =  45  deg,  r/s  ^0.5  results  are  also  used  for  r/s  ^  0.2  except  the  [Kb(W)]min 
lower  than  that  for  «I>  =  0  deg  (see  Figure  24).  The  minimum  values  of  Kb(w)  for  angles  of  attack 
above  the  upper  boundary  shown  in  Figure  26,  where  fm  choking  is  important,  are  twice  this 
minimum  value  for  «1>  =  45  deg  shown  in  Figure  24.  Even  though  this  value  is  twice  the  ^  =  45  deg 
results,  it  is  still  much  less  than  the  ^>  =  0  deg  results.  This  methodology  thus  says  that  the  body¬ 
wing  interference  for  <I)  =  45  deg  at  high  angle  of  attack  is  much  less  than  the  <I>  =  0  deg  results,  but 
when  fin  choking  is  present,  the  body-wing  interference  is  closer  to  the  ^  =  0  deg  results.  In  other 
words,  this  is  an  approximate  method  to  account  for  the  fact  that  the  pressures  in  the  vicinity  of  the 
wing-body  region  are  increased  substantially  when  sonic  flow  is  present.  The  data  for  the  r/s  =  0.5 
results  of  course  come  from  Tables  1-11.  Between  the  two  boundaries  on  Figure  26,  linear 
interpolation  is  used. 

The  other  physical  phenomenon  present  in  the  Figure  26  results  is  an  apparent  "fin-fin 
interference"  effect  that  occurs  at  0  =  45  deg  and  high  AOA  that  is  not  present  for  «!>  =  0  deg.  In 
examining  Table  6,  it  is  seen  that  is  less  than  1 .0  for  low  Mach  numbers.  This  is  based 

on  the  Figure  3  data  and  substantiated  with  validation  against  other  results.  Unfortunately, 
equivalent  high  AOA  missile  component  data  were  not  available  for  supersonic  configurations.  As 
a  result  most  of  the  values  for  [K„_,]  .  were  assumed  to  be  1.0. 

Recall  that  the  way  Kw(b)  was  calculated  was  based  on  the  wing-body  data  for  O  =  45  deg 
and  the  wing  alone  data  for  <I>  =  0  deg.  Hence,  it  is  suspected  that  the  reason  the  value  of 

1.0  at  low  Mach  number  is  really  a  fin-fin  interference  effect  for  fins 
oriented  at  ^  =  45  versus  0  deg  roll.  It  is  suspected  that  the  leeward  plane  fins  are  completely  in 
the  wake  of  a  low  dynamic  pressure  region  (M  =  0. 1)  and  hence  do  not  have  the  normal  force  that 
the  cruciform  fins  have  at  <I>  =  0  deg.  In  other  words,  it  is  suspected  that  a  value  of  less  than  one 
would  be  obtained  for  fins  at  <I>  =  45  deg  versus  fins  at  «I>  =  0  deg,  even  if  no  body  were  present. 
Since  fin-fin  interference  is  not  treated  separately,  it  is  included  in  the  Kw(b)  term. 

The  interference  factor  methodology  for  r/s  ^  0.2  is  thus  designed  to  be  independent  of  roll 
for  small  angles  of  attack  and  to  start  showing  the  roll  dependence  when  "fin  choking"  or  fin-fin 
interference  effects  become  important. 

There  is  also  a  change  in  the  center  of  pressure  calculation  for  low  r/s  value  configurations, 
which  will  be  discussed  in  section  2.2.6. 

0.2  <  r/s  <0.5.  For  r/s  values  between  0.5  and  0.2,  the  interference  factors  (Kw(b)  and  Kb(W)) 
are  computed  for  r/s  values  of  0.2  and  0.5  using  the  approaches  just  outlined,  and  then  the  values 
at  the  given  r/s  of  interest  are  computed  based  on  linear  interpolation.  That  is 
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r/s-0.2 

0.3 


{  [^B(W)L=o.5 


(17B) 


(8  -  lOr/s) 


(17C) 


Here  Kg  is  the  minimum  value  shown  in  Figure  24  represented  by  the  lower  curve  for  the  results 
when  ^>  =  45  deg. 


This  methodology  inherently  assumes  that  the  "fin  choking"  and  low  r/s  phenomena  occur 
in  a  linear  fashion  as  r/s  gets  smaller.  Realistically,  this  is  probably  not  the  case;  however,  only 
numerical  computations  or  wind  tunnel  results  can  specify  the  detailed  fashion  of  how  the 
interference  factors  vary  with  r/s.  Until  such  results  are  available,  or  until  computations  are  made 
on  configurations  outside  the  data  base  upon  which  this  methodology  is  based,  showing  it  to  be 
incorrect.  Equation  (17)  will  be  used. 


2.2.6  Center  of  Pressure  of  Wing  and  Interference  Lift  Due  to  AO  A  at  <I>  =  45  deg 

There  are  two  changes  in  the  center  of  pressure  prediction  for  the  interference  lift  from 
<E>  =  0  deg  to  ^>  =  45  deg  roll  orientation.  Before  the  changes  are  discussed,  a  brief  review  of  the 
=  0  deg  roll  orientation  method  is  appropriate. 

The  AOA  prediction  of  center  of  pressure  of  the  Kw(b)  term  is  based  on  linear  theory  of  the 
wing-alone  lift  at  small  AOA  and  goes  to  the  centroid  of  the  wing  planform  area  at  large  AOA 
(«w  =  60  deg).  This  transition  is  accomplished  in  a  quadratic  fashion.  If  A  and  B  are  the  centers 
of  pressure  of  the  linear  and  nonlinear  normal  force  terms  (in  percent  of  mean  geometric  chord),  and 
«w  “  +  6,  then  the  center  of  pressure  of  the  wing-body  or  wing  alone  lift  is 


{^cp)^B  (^cp)^  ^  ^  36  A)  °'w  (A  B)  (18) 

Equation  (18)  is  the  methodology  used  for  roll  position  of  0  deg. 

When  the  fins  are  rolled  to  a  non-zero  roll  orientation,  the  center  of  pressure  Equation  (18) 
will  change  because  of  the  geometry  of  the  wings  and  an  asymmetric  effect  on  the  wing  loading. 
To  visualize  this  effect,  imagine  a  missile  rolled  to  ^  =  45  deg  and  increasing  in  AOA.  As  AOA 
increases,  two  things  occur.  First,  the  windward  plane  fins  carry  more  and  more  of  the  load 
compared  to  the  leeward  plane  fins.  Second,  the  local  Mach  number  in  the  windward  plane  is 
different  and  typically  lower  than  the  leeward  plane.  This  has  the  effect  of  .shifting  the  wing  alone 
center  of  pressure  forward  in  the  windward  plane.  Since  the  load  and  wing  centers  of  pressure  are 
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different  on  the  windward  and  leeward  plane  fins,  this  results  in  a  net  forward  shift  in  the  center  of 
pressure  for  <1>  =  45  deg  roll  compared  to  the  ^>  =  0  deg  computation  of  Equation  (18).  This  shift 
appears  to  occur  for  all  Mach  numbers,  and  is  largest  at  moderate  AOA,  and  goes  to  zero  at 
AOA  0  and  90  deg.  At  90  deg,  the  windward  plane  fins  carry  almost  all  the  load  compared  to  the 
leeward  plane  fins,  but  geometrically,  the  fins  are  all  aligned  perpendicular  to  the  AOA  plane. 

Mathematically,  this  geometrical  shift  can  be  approximated  by: 


^  Ycp  cos  O  sm(2«) 


AC 


N, 


W(B) 


'N, 


'W(B) 


(19) 


Ycp  of  Equation  (19)  is  the  centroid  of  the  planform  area  in  the  Y  direction  of  the  wing 
planform.  If  the  wing  is  assumed  to  be  represented  by  a  trapezoidal  shape,  then 


Y 


CP 


=  r 


b  \ 

C.  _  c, 

C,  -  c,j 

.2  3) 

(20) 


The  cos  <I>  factor  rotates  the  ACn  term  from  being  normal  to  the  wing  planform,  to  the  plane  of  the 
normal  force  vector.  The  sin  2«  term  allows  the  COP  shift  to  vary  between  0  at  AOA  0  and  90  deg 
to  a  maximum  at  AOA  of  45  deg.  AC.,  /C.,  represents  the  difference  between  the  load  of  the 

^  "^WCB)  ^W(B) 

windward  and  leeward  plane  fins  as  a  fraction  of  the  total  wing-body  load.  Using  the  NASA  data 
base*,  an  estimate  was  made  for  several  wings.  An  approximate  linear  trend  with  angle-of-attack 
was  arrived  at.  While  data  are  not  available  above  «  of  40  deg,  it  will  be  assumed  this  trend  is  linear 
between  AOA  of  0  and  65  deg.  At  «  ^  65  deg  it  is  assumed  the  windward  plane  fins  carry 
80  percent  of  the  load  of  the  wings  and  the  leeward  plane  fins  carry  only  20  percent.  With  this 
assumption. 


AC 


N, 


'WCB) 


"N, 


'W(B) 


«  <  65° 


-  O.g  ;  «  >  65° 


(21) 


In  Equation  (21),  « 


is  in  degrees.  Combining  Equations  (19),  (20),  (21),  there  is 


r  + 


c.  -  c, 


22 

3; 


cos(O)  sm(2«) 


0.8p< 
65  , 


obtained 


«  ^  65 
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=  -  0.8 


r 


C  •+  C 


ty 


cosO  sm(2«)  ;  «  >  65° 


(22) 


Equation  (22)  is  added  to  Equation  (18)  for  the  roll  orientation  of  45  deg. 

The  other  change  being  made  in  the  center  of  pressure  for  =  45  deg  is  to  use  the  local 
Mach  number  as  AOA  is  increased  (versus  the  freestream  Mach  number)  to  compute  the  body 
carryover  COP.  The  local  Mach  number  can  be  computed  by  Equation  (15)  using  MNT.  The  local 
Mach  angle  is  then 


(23) 


Using  the  local  Mach  angle  to  determine  the  angle  over  which  linear  theory  impacts  the  body 
carryover  lift  has  the  effect  of  shifting  this  component  forward.  That  is  because  Ml  decreases  with 
increasing  AOA  causing  //l  to  approach  90  deg.  At  90  deg  Mach  angle,  the  center  of  pressure  of 
linear  theory  is  at  the  same  location  of  slender  body  theory,  that  is,  at  the  centroid  of  the  wing  root 
chord. 


2.2.7  Wing-Bodv  and  Bodv-Wing  Interference  Due  to  Control  Deflection 

The  same  general  approach,  with  slight  modifications,  to  the  nonlinear  model  of  Reference  1 
for  4>  =  0  deg  roll  will  be  used  here  for  the  O  =  45  deg  roll  position.  In  the  Reference  1  method, 
kw(B)  and  kB(W)  were  approximated  by 

^w(B)  ~  Cj(M)  [1^w(B)^bj  ’’’  ^)  (24) 


k 


B(W) 


(25) 


The  parameters  Ci  and  C2  were  derived  based  on  numerical  experiments  of  the  AP95  compared  to 
data.  This  was  because  many  of  the  fins  in  the  Reference  8  data  base  were  too  small  to  allow 
accurate  estimates  of  kw(B)  as  a  function  of  parameters  of  interest.  As  a  result,  total  missile  load  data 
were  used  and  empirical  values  of  Ci  and  C2  estimated  as  a  function  of  combined  local  angle  of 
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attack  of  the  wing  H  +  6|,  and  Mach  number.  Figure  27  gives  the  values  of  the  parameters  Ci  and 
C2  for  =  0,  and  Figure  28  gives  the  complementary  set  for  =  45  deg. 

In  examining  the  constants  and  model  of  Figure  27,  several  physical  phenomena  occur  that 
are  modeled  in  a  semiempirical  sense  by  Equations  (24)  and  (25).  These  phenomena  are 
qualitatively  shown  in  Figure  29 A.  At  low  Mach  number,  Figure  29A  indicates  the  SET  gives  a  low 
value  of  kw(B)  for  small  values  of  .  At  a  value  of  of  about  25  deg,  the  controls  lose 
effectiveness  as  a  result  of  a  combination  of  stall  and  blow-by  effects  due  to  the  separation  between 
the  wing  and  body.  At  an  of  about  55  deg,  the  controls  have  lost  all  effectiveness.  At  Mach 
numbers  greater  than  about  4,  the  controls  initially  generate  less  effectiveness  than  is  generated  by 
SET  for  values  of  «w  up  to  about  20  or  25  deg.  The  controls  then  become  more  effective  because 
of  nonlinear  compressibility  effects.  On  the  other  hand,  at  an  of  around  45  to  50  deg,  the 
controls  once  again  begin  to  lose  effectiveness,  presumably  because  of  shock  interactions  and  blow- 
by  effects.  For  Mach  numbers  in  between  subsonic  and  high  supersonic,  kw(B)  has  behavior  in 
between  the  two  extremes  illustrated  in  Figure  29 A. 

In  comparing  the  nonlinear  control  deflection  models  for  $  =  0  and  45  deg  roll  in  Figures  27 
and  28,  a  lot  of  similarity  is  seen.  The  constants  for  the  ^  =  45  deg  are  slightly  different  than  those 
for  $  =  0  deg  and  the  values  of  where  the  nonlinearities  begin  are  somewhat  different. 
However,  by  and  large.  Equation  (24)  holds  for  both  the  $  =  0  and  45  deg  roll  cases.  It  should  be 
pointed  out  that  in  Reference  1,  mostly  linear  variations  of  kw(B)  with  were  used.  However,  these 
have  been  improved  upon  for  the  $  =  45  deg  case  with  cubic  fits  of  control  deflection  data  as  seen 
in  Figure  28.  As  such,  all  nonlinear  effects  are  included  in  the  variations  of  kw(B)  as  a  function  of 
Mach  number  and  | «  +  6 1 . 

Figure  29E  assumes  that  kB(W)  can  be  represented  by  SET  up  to  some  value  of  |  «vv  |  >  at  which 
point  it  decays  to  a  percent  of  SET.  For  $  =  0  deg,  these  values  are  |  I  =70  deg  and  50  percent 
of  SET  as  a  minimum.  For  $  =  45  deg  roll,  the  model  for  kB(w)  is  based  on  |6|  only  and  begins 
decaying  at  jS  |  =  0.  It  reaches  a  minimum  of  kB(w)  of  25  percent  of  SET  analogous  to  the  [KB(W)]min 
of  Figure  24. 

It  should  also  be  noted  that  kw(B)  and  kB(w)  of  Figure  28  are  multiplied  by  1.414  to  indicate 
that  all  four  fins  are  assumed  to  be  deflected  by  an  equal  amount  in  the  $  =  45  deg  roll  position. 
Finally,  for  Mach  numbers  in  between  the  values  on  Figures  27  and  28,  linear  interpolation  is  used. 


2.2.8  Nonlinear  Wing-Tail  Interference  Model  for  O  -  0  deg 

The  nonlinear  wing-tail  interference  model  used  for  the  $  =  0  roll  is  defined  by 


(^NocIt 


^  ^W(B>  ^  (®T  ’’t)  -^W 
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TCV) 


2^(AR)t  (^w  ’•w) 
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MACH 


NONLINEAR  MODEL 


M^O.8 


M=l.l 


If  I  “tv  I  ^  24.0-‘kTv(B)  ~  Ll[k,V(B)]sB 

If  j“Tv|  >24.0-kw(B)  =  L1[0.000794|«tvP  -  0.0933  |“tv|  +2.71] 

F  =  1.I 


If  1  “tv I  ^  15.0-*kTv(B)  ~  L0|kTv(B)3sB 

If  i“Tvi  >  15.0-kTv(B)=  1.0[0.00087|“tv|'  -  0.0825|“tv|  +  1.98] 
F  =  l.l 


M  —  1.5  If  |«Tv|  ^  l0.0->kTv(B)  —  LOPctv(B)]sb 

If  i“wi  ^  10.0-‘kTV(B)  =  L0Pctv(b)]sb  ~  0.005[|“tv|  -  10.0] 

If  |“tvU  20.0-F  =  0.8 

If  I  “tvI  >  20.0-*F  =  0.8  +  0.10[|  «Tv|  -  20.0] 


M— 2.0  If  |“tv|  ^  10.0-*kTv(B)  — 0.9|kTv(B)]sB 

ff  |‘*wi  ^  10.0-*kTv(B)  ~  0.9[kTv(B)]sB  ~  0.003[|“tv|  “  10.0] 

If  j“Tvi  ^20.0-F  =  0.8 

If  i“wi  >20.0-F  =  0.8  +  0.17[|cc^|  -  20.0] 


If  |“tv|  ^  40.0-*kTv(B)  “  0.9PCtv(b)]sB 

ff  |“w|  ^  40.0->kTV(B)  ~ 0.9[kTv(B)]sB  +  0  005[| «Tv|  ~  40.0] 

If  i“Tvi  s:  30.0-F  =  0.9 

If  |“Tvi  >30.0-F  =  0.9  +  0.15[|<»:^|  -  30.0] 


M  —  2.87  If  |“tv|  ^  40.0-‘k'^(B)  —  0.9[k^(BjsB 

^i“wi  ^  40.0-*kTV(B)  ~ 0.9jkTV(B)]sB  +  0.005[| “tvI  ~  40.0] 

If  j«Tvi  ^  30.0-F  =  0.9 

If  I “tvI  >30.0-F  =  0.9  +  0.1711  “tvI  "  30.0] 


M  — 3.95  If|“vi;|  ^  20.0-*k^(B)  —  0.8(k^(B)]5B 

If  I  “w|  ^  20.0-*kTV(B)  ~  0.8(kTV(B)]sB  +  0.007[|  “tvI  ~  20.0] 

If  |“tv|  ^  30.0-F  =  0.9 

If  |«^j  >30.0-F  =  0.9  +  0.2[|“,v|  -  30.0] 


M  ^  4.6  If  |“tv|  ^  20.0-*kTV(B)  “  0.75[kTV(B)]sB 

If  i“wi  ^  20.0-*kTV(B)  =  0.75|kTv(B)]sB  +  0.013[| “tv|  “  20.0] 

If  l“,vi  ^  35.0-F  =  0.9 

If  jo^Tvi  >35.0-F  =  0.9+0.2[j“Tv|  -  35.0] 


Where  o^tv  =  “  +  ^  and  for  large  «tv 


l^Bcw)  ~  [I^B(w)^sB  for  |“w|  ^70  deg 


*^BCW)  II^B(W)IsB  I 


l^CB)  I*^cb)IhwI  =70 


for  70  <  |“wl  s  90 


“  -70  1 

LJL -  for70<|-Tvl  ^90 


I^Bov)  “  IrTV(B)  ~  0  for  I  “tv  I  >  90  deg 


FIGURE  27.  SEMIEMPIRICAL  NONLINEAR  CONTROL 
DEFLECTION  MODEL  ($  =  0  DEG) 
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NONLINEAR  MODEL 


If  |“w|  <.  40.0  deg-*kTj,(B)  “ 

If  >40.0  deg-*kw(B)  =  L15[0.000794|«wP-  0.0933  |«w|  +2.71] 
F=L1 


<.  15.0-‘k^(B)  —  0.95[k^(B)]sB 

>  15.0-kw(B)  =  0.95[0.00087|«w|'  -  0.0825 |«w|  +  1-98] 


If  I  “^1  i  35.0-*k^(B)  —  0.95[k^(B)]sB 

If35.0  <  l-wl  ^  55.0-kw(B,=  [8.067  x  lO’’  (|-^|  -  35.0)^  +  0.00201  (|-w|  -  35.0)^ 

-0.0295  (l-w!  -  35.0)  +  0.939] 

hml  -  55.0 

If55.0<lc-w|  ^  62.0 -*kw(B)=  1.02  1.0  -  -  [kw(B)]sB 

7  v2S4 

^f  |“wl  ^  62.0-*k^(B)  =  0.04[kw(B)]sB 
If  j«wi  ^  20.0-F  =  0.8 

If  >20.0-F  =  0.8  +  0.17[l«wl  -  20.0] 


If  |“^|  ^  20.0-*k^(B)  ~  0.95[k^(B)]sB 

If  20.0  <  l-xwl  ^  55.0-kw(B)=  [5.699  x  lO'*  (l«w|  -  20.0)^  +  0.00420  (|«w|  -  30.0)^ 

+9.00583  (|«w|  -  20.0) +  0.980]  [kw(B)]sB 

l“n,l  -  55.0 

If55.0<|«wU62.0-kw(B)  =  0.89  1.0  -  -  [kw(B)]sB 

7.3279 

If  I  “wi  ^  62.0-*kTj,(B)  “  0.04[kix,(B)]sB 

If  ^  20.0-F  =  0.8 

If  |«^|>20.0-F  =  0.8  +  0.17[KI  -20.0] 


If  I'Xwl  ^  30.0-*k^(B)  -  0.95|k^(B)]sB 

If  30.0  <  |«wh  55.0-kw(B)=  [4.257  x  lO*’  (|«^|  -  30.0)’  +  0.00291  (|«w|  "  30.0)’ 

+0.0388  (|«^|  -  30.0)  +  0.976]  \k^(s)U 

|«  I  -  55.0^ 

If  55.0  <  I  «w|  ^  90.0  -kw(B)  =  [kv(B)].=55  -  -  ([kw(B)]>=55 

V  / 

-  0.04[kw(B)]sB) 

If  I«w|  ^30.0-F  =  0.9 

If  I >  30.0-F  =  0.9  +  0.15[|«^|  -  30.0] 


If  [“‘wl  ^  30.0-*k^(B)  —  0.95[kw(B)]sB 

If  30.0  <  |«w|  ^  55.0-kw(B)=  [6.526  x  10  ’  (j«^|  -  30.0)’  -  0.00405  -  30.0)’ 

+9.0575  (|«w|  -  30.0)  +  0.947]  [kw(B)]sB 

(  |«  I  -  55.0 

If  55.0  <  |“wl  ^  90.0-*k^(B)=  [k^fBvlo,=s5  “  - rrr -  ([kw(B)]<=55 


If  55.0  <  |“w|  ^  90.0 ->k^(B)  [IS(r(B)]o<=55  ' 

-  0.04[k^(B)]sB) 

If  I«wi  ^  30.0-F  =  0.9 

If  >30.0-F  =  0.9  +  0.17[|«w|  -30.0] 
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MACH 

NONLINEAR  MODEL 

M  =  3.95 

If  i“w|  ^  —  0.88[kTiy(B)]sB 

If  35.0  ^  |«w|  ^  55.0-kw(B)=  [-8.84  x  IQ-'  ( 

+0.0397  (l-wl 

If  55.0  i  |«^|  i  90.0-*k^(B)= 

-  0.04[kTj,(B)]sB 

If  ^  30.0-F  =  0.9 

If  >30.0-F  =  0.9  +  0.2[|c.^|  -  30.0] 

«w|  -  35.0)^  +  0.000173  (|<=^tj,|  -  35.0)^ 

-  35.0)  +  0.884]  [kTj,(B)]sB 

-  55.o) 

[  35.0  j 

) 

Ms:  4.6 

If  ^  35.0-*k^(B)  =  0.83[k^(BdsB 

If  35.0  s  |«w|  ^  55.0-kw(B)=  [4.697  x  l0-’( 

+0.0739  (l-wl 

If  55.0  ^  1“^!  <,  90.0-»k^(B)“  [^v(B)l«=ss 

-  0.04[k^(BdsB 

If  l-wl  ^35.0-F  =  0.9 

If  >35.0-F  =  0.9  +  0.2[|ccw|  -35.0] 

|«w|  -35.0)’- 
-  35.0) +  0.83 

'  l“wl  -  55.0 1 

1  35.0  J 

) 

0.00463  (1^,^1  -35.0)’ 

9]  nrw(B)]sB 

(P^W(B)]«=5S 

0  ^  M  ^  » 

|Sn,l 

kB(w)  “[ksonolsB  ^  (0.75  [kflfvsolsB)  f®r  1^1  ^  30.0 

kBov)  “  0.25  [kB(w)]sB  for  |6w|  ^  30.0 

FIGURE  28.  SEMIEMPIRICAL  NONLINEAR  CONTROL 
DEFLECTION  MODEL  (4&  =  45  DEG)  (CONTINUED) 


where  kw(B)  and  F  were  defined  in  Figure  27.  Also  Cj^  of  both  the  wing  and  tail  were  nonlinear 
values.  The  interference  factor  i  was  defined  by  the  traditional  SBT/LT^*  approach. 

While  Equation  (26)  was  helpful  in  defining  some  of  the  nonlinearities  of  missile 
aerodynamics  not  predicted  by  SBT  or  LT,  it  still  had  two  major  weaknesses.  First  of  all,  the  tail 
interference  factor  methodology  was  independent  of  Mach  number,  and  second,  there  was  no 
parameter  for  the  first  term  of  Equation  (26),  similar  to  the  parameter  F  of  the  second  term,  to  help 
control  the  quasilinear  character  of  the  equation.  These  two  problems  were  partially  accounted  for 
in  the  AP95*  by  placing  upper  limits  on  theC„  term  as  a  percent  of  total  tail  lift  and  as  a  function 

^TCV) 

of  Mach  number.  However,  these  upper  limits  mainly  affected  large  wing  configurations  where  the 
wing-tail  interference  could  exceed  the  tail  lift. 

In  an  effort  to  improve  upon  the  present  semiempirical  wing-tail  methodology,  the  literature  , 
was  searched  for  wing-tail  interference  data.  References  22  and  23  were  found  in  this  process. 
These  references  documented  wind  tunnel  tests  on  two  different  missile  configurations  in  which 
actual  wing-tail  normal  force  measurements  were  made.  Reference  22  measured  these  results  on 
individual  fins  as  a  function  of  roll  position  at  Mach  numbers  of  2,  3  and  4.  Reference  23  gave 
results  for  only  the  $  =  0  deg  roll  position  at  M„  =1.1.  Figure  30  shows  the  configuration  tested 
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CONTROLS  BECOME  MORE 
EFFECTIVE  DUE  TO  COMPRESSIBILITY 
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FIGURE  29A.  QUALITATIVE  TREND  OF  WING-BODY  INTERFERENCE  DUE 
TO  CONTROL  DEFLECTION  AS  FUNCTION  OF  M., 
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FIGURE  29B.  QUALITATIVE  TREND  OF  BODY-WING  INTERFERENCE  DUE 
TO  CONTROL  DEFLECTION  AS  FUNCTION  OF  cc^ 
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MS  45.705 


A.  CONFIGURATION  TESTED  IN  REFERENCE  22 


B.  SINGLE  FIN  WING-TAIL  INTERFERENCE  (M=L96) 


C.  SINGLE  FIN  WING-TAIL  INTERFENCENCE  (M=4.02) 


FIGURE  30.  CONFIGURATION  AND  SINGLE  FIN  DATA  FROM  REFERENCE  22 
FOR  WING-TAIL  INTERFERENCE  NORMAL  FORCE 
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in  Reference  22  along  with  results  at  =  1.96  and  4.02  at  roll  position  of  0  deg.  Also  reported 
in  Reference  22  were  the  theoretical  computations  from  analytical  codes  referred  to  as  AERODSN,^'* 
MISSILE  2A,^^  and  LRCDM,^®  According  to  Reference  22,  the  AERODSN  code  has  the  same 
wing-tail  interference  model  as  the  older  version  of  the  Aeroprediction  Code  and  is  based  on  SET 
froin  Reference  21.  The  MISSILE  2A  program  is  based  on  theory  and  experiment  and  therefore 
should  contain  some  of  the  nonlinearity  associated  with  the  wing-tail  interference.  The  LRCDM 
code  is  based  on  paneling  methods  and  data  bases  and  therefore  should  also  contain  some  of  the 
nonlinearities  of  the  wing-tail  interference.  For  comparison  purposes,  computed  by  the  AP95 

^T(V)  ^ 

for  a  single  fin  is  also  shown  on  Figures  30B  and  30C.  The  AP95,  although  it  contains  some 
nonlinearities,  still  resembles  the  SET,  as  seen  from  the  close  proximity  to  the  AERODSN  results. 
It  is  fair  to  say,  in  viewing  the  comparison  of  the  various  theoretical  approaches  for  predicting  wing- 
tail  interference  in  Figure  29,  that  improvements  in  the  theory  are  needed. 

One  might  argue  that  the  value  of  is  so  small  that  it  can  be  neglected.  However,  in 
most  cases  the  moment  arm  is  fairly  large  so  even  a  small  loss  of  lift  on  the  tail  is  magnified  in  the 
pitching  moment.  As  a  result  of  this  and  the  fact  the  Cj^  term  is  not  predicted  well  by  any  of  the 
theoretical  approaches,  it  was  decided  to  try  to  develop  a  new  semiempirical  approach  to  predict 
.  The  new  approach  will  still  be  anchored  by  the  traditional  SET^*  but  will  be  modified 

Mt(v) 

according  to  the  data  from  References  22  and  23. 

A  model  for  no  control  deflection,  based  on  a  third  order  equation  in  AOA  is  postulated  to 
fit  the  data  from  both  References  22  and  23.  This  model  is  defined  by 

=  A  +  (27) 


The  four  constants  of  Equation  (27)  require  four  conditions  to  define  them.  These  conditions  are 
as  follows: 

(1)  Slope  of  SET  value  for  Cjj  near  «  =  0  (Using  this  value  and  a  value  of  the  real 
slope  from  data,  a  modifieTvalue  of  this  slope  can  be  obtained.) 

(2)  AOA  where  Cj^  goes  to  zero  based  on  data 

(3)  AOA  where  C„  reaches  a  maximum  based  on  data 

^TCVL 

(4)  Maximum  value  of  as  a  percent  of  SET  at  AOA  whereC^^  is  a  maximum 

Note  in  the  four  conditions  chosen,  SET  is  used  twice.  Eased  on  these  conditions  the  four  constants 
of  Equation  (27)  become 


A  =  0 


B  = 


Nt(v) 

1  J 

o<=0 

SBT 
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C 


-B-  D 


oc 


2 

N 


oc 


N 


D 


<  —  Roc  oc  +  B 

N  “  N  F  “ 


3  2  2 

O'  ~  O'— 

N  F  F  N 


oc 


2 

F 


The  parameters  in  the  constants  B,  C,  D 
O'  =  Value  from  Figure  3 1 A 
eop  =  (Value  from  Figure  3  IB)  x 
El  =  Value  from  Figure  3 1C 
El  =  (Value  from  Figure  3  ID)  x 


are  defined  as  follows: 


100 


Referring  back  to  Equation  (26),  if  we  break  this  equation  down  into  the  part  due  to  AOA  and  the 
part  due  to  control  deflection,  we  have 


+ 


(28) 


Equation  (27)  is  used  to  estimate  the  term  due  to  oc,  and  the  section  2.2.4  methodology  is  used  to 
predict  the  methodology  for  the  second  part  of  Equation  (28).  We  now  have  independent  controls 
on  both  portions  of  Equation  (26),  which  is  based  on  SBT,  but  the  controls  allow  the  SBT  to  be 
adjusted  to  data. 

The  question  arises  regarding  how  to  account  for  fins  that  are  different  sizes  and  locations 
than  those  tested  in  References  22  and  23 .  Fortunately,  both  fins  tested  in  these  references  were 
significantly  larger  than  the  tail  surfaces  and  were  located  in  a  wing  versus  a  canard  location. 
Hence,  the  present  approach  will  be  to  use  these  results  directly  for  wings  or  canards  of  less  area  to 
body  reference  area  than  those  tested.  For  wings  of  greater  area,  it  is  assumed  that  the  AOA, 
where  theC„  becomes  negligible  is  increased  according  to 


=  O',.,  for  A„,/A„„„  <  5.5 


OC  =  OC 

N 


^  A  /A  ^ 

REF 


5.5 


for  A^/A^p  >5.5 


(29) 


with  an  upper  limit  on  of  2. 5  “'jj  . 
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FIGURE  3 1  A.  AOA  WHERE  WING-TAIL  INTERFERENCE  IS  NEGLIGIBLE 


FIGURE  3  IB.  AOA  WHERE  WING-TAIL  INTERFERENCE 
IS  A  MAXIMUM  (PERCENT  OF 
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FIGURE  3  ID.  SLENDER  BODY  THEORY  PREDICTION  OF  WING-TAIL 

INTERFERENCE  AT  AOA  WHERE  [C„  1  REACHES  A 

tCV)JexP 

MAXIMUM  AS  A  FRACTION  OF  EXPERIMENTAL  DATA 
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The  value  of  A^/Areij  of  5.5  corresponds  to  the  wing  area  ratio  of  Reference  22.  Also,  for 
wings  larger  than  those  of  Reference  22,  an  upper  limit  on  the  amount  of  lift  loss  on  the  tail  will 
remain  in  effect.  This  upper  limit  is  defined  by  the  following  methodology: 


=  1.0  ;  «<5 

=  1.0  -  0.04125  («-5)  ;  «  >  5  (29A) 


=  0.9  -  0.025«  ;  «  <  10 

=  0.65  -  0.0235  («  -  10)  ;  «  >  10  (29B) 


=  0.8  -  0.025«  (29C) 

where  «  is  AO  A  in  degrees.  Equation  (29  A)  says  that  at  «  =  0  deg,  the  maximum  lift  loss  on  the 
tail  is  limited  to  100  percent  of  the  tail  lift,  regardless  of  the  size  of  the  wings.  The  percent  lift  loss 
then  decreases  linearly  with  AO  A  as  defined  by  Equations  (29  A)  -  (29C). 

Admittedly,  this  is  conservative  (overpredicts  )  for  values  of  A^/Aref  <5.5  and  is 
simply  a  judgement  based  on  numerical  experiments  for  values  of  Aw/Aref  >5.5.  It  does 
accomplish  the  objective  of  making  the  wing-tail  interference  with  no  control  deflection  more 
closely  approximate  data  than  available  approaches,  including  the  AP95.  This  is  illustrated  by  the 
results  shown  in  Figures  SOB  and  30C  for  Improved  Aeroprediction  Code  1995  (IAP95). 

Before  moving  to  the  wing-tail  interference  methodology  for  the  $  =  45  deg  roll  position, 
a  comment  would  be  valuable  on  the  Figure  3 1  results,  which  basically  compare  SBT  to  data.  First 
of  all,  it  is  clear  that  at  low  AO  A,  SBT  underpredicts  for  low  Mach  number  and  overpredicts 
it  at  high  Mach  number,  for  the  Reference  22  configuration.  The  point  of  optimum  prediction 
appears  to  be  around  Mach  2  (see  Figure  3 1C).  Second,  the  C„  term  decays  much  faster  at  high 

■^TCV) 

AOA  than  does  SBT.  This  is  increasingly  true  as  Mach  number  increases.  This  again  highlights 
the  Newtonian  Impact  assumptions  at  high  Mach  number  where  any  vortices  in  the  leeward  plane 
are  completely  dominated  by  the  dynamic  pressure  in  the  windward  plane.  Figures  3 1 A  and  3  IB 


ForM^  1.5 


Cm 

^TCV) 


For  1.5  <M^  2.5 


Cm 

^T(V) 


For  M>  2.5 
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illustrate  this  fact,  showing  that  the  AOA  where is  negligible  gets  smaller  as  M„  increases; 

T(V)  , 

also,  the  maximum  magnitude  as  a  percent  of  SET  gets  smaller  with  increasing  M„. 


2.2.9  Nonlinear  Wing-Tail  Interference  Model  for  ^  =  45  deg 

The  SET  representation  of  the  wing-tail  interference  for  $  =  45  deg  roll  is  outlined  in 
Appendix  A.  As  pointed  out  in  the  appendix,  the  assumption  is  made  that  the  strength  of  the  vortex 
shed  from  the  windward  plane  fins  is  equal  in  magnitude  and  opposite  in  direction  to  the  leeward 
plane  fins.  At  small  AOA,  this  assumption  is  quite  reasonable  and  is  partially  what  leads  to  the  fact 
that  cruciform  wing-body-tail  missile  aerodynamics  are  independent  of  roll  position.  As  AOA 
increases,  this  assumption  becomes  less  and  less  valid.  In  fact,  the  lift  of  the  windward  plane  fin 
is  much  larger  than  the  leeward  plane  fm  as  AOA  increases.  This  was  modeled  approximately  in 
the  center  of  pressure  shift  discussed  earlier  by  approximating  a  linear  variation  in  the  shift  of 
normal  force  to  the  windward  plane  fin  from  the  leeward  plane  fin  up  to  AOA  65  deg.  At  that  point, 
the  ratio  of  the  windward  to  leeward  plane  load  remained  constant. 

If  we  define  the  factors, 

P„  =  (l.O  +  0.6^j.  P,  =  |l.O  -  0.6.^j  ;  «  s  65 


P„  =  1.6  .  P,  =  0.4  ;  «  >  65 


(30) 


then  the  interference  factors  can  be  weighted  by  Equation  (30)  depending  on  where  the  vortex  is 
shed.  If  it  is  shed  in  the  windward  plane,  then  the  Pw  factor  is  appropriate;  whereas  if  it  is  shed  from 
the  leeward  plane,  then  the  Pj  factor  in  Equation  (30)  is  appropriate.  This  is  an  approximate  way 
to  represent  the  nonlinear  nature  of  the  load  in  the  windward  to  leeward  plane  fm,  and  the  strength 
of  the  vorticity  shed  from  each  fm. 

Using  Equation  (30)  along  with  Equations  (A-42)  and  (A-43)  of  the  appendix,  an 
approximate  nonlinear  representation  of  the  wing-tail  interference  factors  for  the  $  =  45  deg  roll 
position  is 
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All  the  nomenclature  in  Equations  (31)  and  (32)  are  defined  in  the  Appendix  A.  The  wing-tail 
interference  normal  force  is  then 


■'^w(^No<)w  (^Noc)t  ^ 


K 


W(B) 


'N, 


TCV) 


27t  (AR)t  -  r^)  A. 


REF 


(sT-rx)  [ii  cosO  +  ^33^ 


For  small  AOA  and  no  control  deflection,  Equation  (33)  reverts  back  to  SET. 

The  final  modification  for  inclusion  of  nonlinearities  into  the  wing-tail  interference  model 
for  the  $  =  45  deg  roll  position  is  to  adjust  Equation  (33)  based  on  experimental  data  analogous  to 
the  $  =  0  deg  methodology  of  Figure  31.  For  the  45  deg  roll  position,  the  only  data  found  on 
individual  fins  were  from  Reference  22.  Figure  32  gives  theC„  values  for  M„  values  of  1 .96, 
3.01  and  4.02.  In  developing  an  analogous  model  to  Figure  3 1  for  $  =  45  deg,  qualitative  use  will 
be  made  of  the  $  =  0  deg  results  at  lower  Mach  number  since  at  least  one  set  of  data^^  existed  for 
M„  =  1.1.  These  results  will  be  used  to  compare  trends  of  data  as  a  function  of  Mach  number  and 
AOA,  not  their  absolute  values. 

Note  the  curves  of  Figure  32  show  a  point  of  inflection  in  the  experimental  data  between  10 
and  15  deg  AOA.  This  is  because  at  very  low  AOA,  the  windward  and  leeward  plane  vortices  shed 
from  the  wings  adversely  impact  the  tail  normal  force.  However,  at  a  slightly  higher  AOA,  the 
windward  plane  wing-shed  vortex  has  a  positive  effect  on  the  leeward  plane  tail  because  of  the 
counterclockwise  vortex  hitting  the  windward  side  of  the  leeward  tail  surface.  As  the  AOA  is 
increased  higher,  the  wing-shed  windward  plane  vortex  rises  above  the  leeward  plane  tail  fm,  at 
which  point  both  wing-shed  vortices  again  have  an  adverse  impact  on  the  tail  normal  force. 

While  a  model  such  as  that  derived  for  the  $  =  0  roll  position  could  be  derived  for  the 
$  =  45  deg  plane,  it  is  more  difficult  because  of  a  lack  of  data  below  M„  =  2.0  and  the  shape  of  the 
curves  in  Figure  32.  As  a  result,  modified  SET  was  used  to  calculate  at  various  AO  As  and 

at  the  three  Mach  numbers  where  data  were  available.  The  results  were  compared  to  the 
experimental  data  of  Figure  32,  and  the  semiempirical  model  of  Figure  33  was  defined  for  no 
control  deflection. 

Referring  to  Equation  (28),  the  first  part  of  the  equation  for  $  =  45  deg,  6w  =  0  is 
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Exp  (ref.  22) 
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FIGURE  32.  WING-TAIL  INTERFERENCE  (0  =  45  DEG,  8^  =  0  DEG) 
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FIGURE  33.  WING-TAIL  INTERFERENCE  MODEL  FOR  NO  CONTROL 
DEFLECTION  AT  ij*  =  45  DEG 
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with  the  same  maximum  constraints  placed  on  Equation  (34)  as  for  the  «]>  =  0  deg  roll  methodology. 
This  constraint  was  given  by  Equation  (29).  Equation  (29A)  basically  allows  the  negative  tail 
normal  force  to  not  exceed  a  certain  percent  of  the  normal  force  of  the  tail  regardless  of  SBT 
predictions.  Equations  (34)  and  (29),  in  conjunction  with  Gj  from  Figure  33,  and  modified  SBT 
[Equations  (3 1)  -  (33)]  defines  the  wing-tail  interference  model  used  for  $  =  45  deg. 


3.0  RESULTS  AND  DISCUSSION 


The  new  nonlinear  aeroprediction  methodology  for  the  =  45  deg  roll  position  has  been 
validated  against  many  configurations  within  and  outside  the  data  bases  upon  which  it  was 
developed.  The  primary  data  bases  upon  which  the  methodology  was  based  was  given  by 
References  6-8.  If  data  points  are  taken  at  angles  of  attack  of  10,  15,  20,  25,  30,  35,  and  40  deg;  at 
Mach  numbers  of  0.6,  0.8,  1.2,  1.5,2  0,  2.5,  3.0,  3.5,  and  4.5;  at  aspect  ratios  of0.5(A.=0,  0.5,  1.0) 
and  1 .0;  a  total  of  227  points  are  available  to  average.  This  is  much  lower  than  in  previous  data 
bases  due  to  omission  of  the  AR  =  2.0  and  4.0  data.  These  data  were  omitted  due  to  the  very  small 
wings  used.  Figure  34  compares  the  new  O  =  45  deg  theory  to  the  data  base  in  terms  of  average  Cn 
and  Xcp  accuracy.  The  average  accuracy  is  defined  as: 


1  ^  I  Theory  -  Exp 

227  1  Exp 


(35) 


As  seen  in  Figure  34,  the  average  accuracy  of  the  $  =  45  deg  methodology  compared  to  the 
data  is  somewhat  worse  than  the  $  =  0  deg  methodology.  However,  this  is  somewhat  misleading. 
This  is  because  the  $  =  0  deg  results  were  taken  directly  from  Reference  1  without  recomputing  the 
results  with  the  revised  "drag  bucket"  discussed  in  section  2.2.2.  Since  the  revised  "drag  bucket" 
is  not  optimized  to  the  Reference  8  database,  it  gave  errors  at  the  lower  Mach  numbers  substantially 
larger  than  those  in  Reference  1.  This  was  the  main  reason  for  the  2.2  percent  and  0.003{  higher 
average  error  on  Cn  and  Xcp,  respectively,  in  Figure  34.  Even  with  these  increases  in  average  errors, 
the  overall  average  of  6.2  percent  on  Cn  and  1 .2  percent  of  body  length  for  Xcp  are  well  within  the 
average  error  goals  of  ±10  percent  for  Cn  and  ±  4  percent  of  body  length  for  center  of  pressure.  It 
is  expected  that  comparisons  on  configurations  outside  the  Reference  8  data  base  at  $  =  45  deg  will 
be  as  good  as  or  better  than  those  in  Reference  1  at  $  =  0  deg. 

Configurations  considered  in  the  validation  process  include  those  with  both  one  and  two  sets 
of  lifting  surfaces.  Those  with  only  one  set  are  either  body-wing  or  body-tail  configurations.  Those 
with  two  sets  are  either  canard-body-tail  or  wing-body-tail.  The  configurations  with  two  sets  of 
lifting  surfaces  will  be  useful  in  validating  all  the  new  methodology  including  the  wing-tail 
interference  with  and  without  control  deflections. 
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DATA  BASK  PARAMETERS 

0.6  ^  ^  4.6 

0.25  ^  AR  ^  4 
0  ^  A,  1.0 
r/s  =  0.5 

BODY:  12.33  CALIBER  OGIVE  CYLINDER 
3  CALIBER  TAN  OGIVE  NOSE 

COMPUTATIONAL  RESULTS  FOR  227  BODY-TAIL  CASES 

«  =  10, 15,  20,  25,30,  (35,40) 


=  0  DEG 

4»  =  45  DEG 

<I»  =  0  DEG 

4»  =  45  DEG 

AP95 

IAP95 

AP95 

IAP95 

AVG ERROR 

4.1% 

6.3% 

.009t 

.012f 

MAX  ERROR 

11.8% 

16.0% 

.04  If 

.030f 

FIGURE  34.  COMPARISON  WITH  EXPERIMENTAL  RESULTS  OF 
REFERENCE  8  DATA  BASE 

3.1  WING-BODY  OR  BODY-TAIL  CONFIGURATIONS 

Seven  configurations  with  one  set  of  lifting  surfaces  are  considered  to  validate  the  new 
methodology  over  a  range  of  angles  of  attack,  Mach  numbers  and  at  both  the  <>  =  0  and  45  deg  roll 
orientation.  Some  of  the  =  0  deg  results  were  given  previously  in  Reference  1,  but  are  shown  in 
conjunction  with  the  O  =  45  deg  results  for  comparison  purposes.  This  allows  the  reader  to  view 
the  theoretical  comparisons  of  the  roll  orientations  simultaneously,  as  well  as  the  impact  of  roll  on 
the  static  aerodynamics.  Emphasis  in  all  the  validations  will  be  on  normal  force,  pitching  moment 
and  center  of  pressure,  as  it  is  well  known  that  the  axial  force  accuracy  of  all  the  various  versions 
of  the  NSWC  Aeroprediction  Code  have  been  very  good.  Some  axial  force  comparisons  will  be 
shown  later  for  configurations  with  two  sets  of  lifting  surfaces. 

The  first  configuration  is  the  wing  B  configuration  of  Figure  2C.  This  configuration  is  the 
one  used  to  determine  the  Kw(B)  and  nonlinear  interference  factors  at  low  Mach  number. 

Normal  force  coefficient  and  center  of  pressure  computations  compared  to  data^  are  presented  in 
Figure  35.  Reasonable  accuracy  levels  are  obtained  for  this  overall  configuration.  If  one  averages 
the  normal  force  and  center  of  pressure  errors  at  each  5  deg  AOA  according  to  Equation  (35), 
average  accuracy  levels  of  10.2  percent  on  and  1.4  percent  of  body  length  for  X(^-p  are  obtained. 
This  accuracy  level  is  at  the  high  end  of  the  ±10  percent  goal  on  normal  force,  but  well  within  the 
±4  percent  I!  on  X^p. 
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FIGURE  35.  NORMAL  FORCE  COEFFICIENT  AND  CENTER  OF  PRESSURE 
COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(CONFIGURATION  OF  FIGURE  2C;  M«  =  O.I) 
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The  accuracy  on  normal  force  was  intentional  at  higher  AO  A.  Note  in  Figure  35  that  the 
predictions  are  consistently  lower  than  the  data  above  AOA  approximately  25  to  30  deg.  This  is 
because  this  configuration  was  tested  with  a  base  mounted  sting.  Reviewing  some  of  the  wind 
tunnel  model  support  interference  literature  (References  27-30),  several  conclusions  were  reached. 
These  were  that:  a)  For  low  Mach  number,  the  problem  of  estimating  interference  effects  of  a  strut 
or  sting  mount  on  the  model  aerodynamics  at  high  AOA  are  not  known  precisely.  (This  is  still  true 
today),  b)  The  preferable  mount  between  a  sting  and  a  strut  at  high  AOA  is  the  sting,  c)  The  sting 
tends  to  give  positive  interference  (C^  too  high)  and  the  strut  negative  interference  (Cn  too  low). 

The  affected  region  seems  to  be  in  the  AOA  range  30  to  80  deg.  Sting  values  can  be  high 
as  much  as  10-15  percent  and  strut  Cn  values  for  struts  mounted  in  the  mid-body  region  can  be  low 
by  as  much  as  25-30  percent.  For  these  reasons,  the  IAP95  methodology  is  intentionally  designed 
to  underpredict  normal  force  on  configurations  at  subsonic  Mach  numbers  at  high  AOA  where  test 
data  are  from  a  sting  mount,  and  therefore  it  is  also  expected  the  predictions  will  be  higher  than 
those  from  a  strut-mounted  model.  The  under/over  prediction  problem  on  Cn  at  high  appears  to 
go  away  at  Mach  numbers  slightly  greater  than  one.  This  is  suspected  to  be  due  to  the  reduction  in 
the  upstream  influence  of  the  sting  on  the  body  and  the  reduction  of  the  wake  effect  on  the  body 
(strut  mount). 

The  second  configuration  considered  is  a  body-tail  configuration  (Figure  3  6 A),  also  at  a  low 
Mach  number  of  0. 15.  Test  results  for  this  configuration  were  reported  in  Reference  3 1,  but  no 
mention  was  made  of  how  the  configuration  was  mounted  in  the  tunnel  (sting  or  strut).  The  test 
results  are  given  in  Figure  37.  It  is  interesting  to  note  that  the  experimental  data  for  Cn  reaches  a 
maximum  around  an  AOA  of  55  deg  and  declines  slightly  up  to  90  deg  AOA.  This  is  the  typical 
behavior  for  a  strut  mounted  configuration.  The  average  errors  for  Cn  and  X^p,  between  the 
predicted  and  experimental  values,  are  10.1  percent  and  4.4  percent  J,  respectively.  Both  of  these 
are  at  the  high  end  of  the  AP  performance  goal  of  ±10  percent  on  Cn  and  ±4  percent  f  on  Xcp.  It  is 
interesting  to  note  that  the  Cn  error  in  this  case  is  positive,  whereas  for  the  case  of  Figure  3  5  it  was 
negative.  These  two  figures  (Figures  35  and  36)  illustrate  the  difficulty  in  accurately  predicting 
static  aerodynamics  at  subsonic  Mach  numbers  for  AOA  greater  than  roughly  30  deg,  because  it  is 
not  clear  what  is  prediction  error  and  what  is  wind  tunnel  measurement  error. 

The  third  body-tail  case  even  more  vividly  illustrates  the  problem  of  predicting 
aerodynamics  at  subsonic  and  transonic  Mach  numbers  at  high  AOA.  The  configuration  is  shown 
in  Figure  36B  and  was  tested  at  Arnold  Engineering  and  Development  Center  (AEDC)  with  a  strut 
mounted  to  the  leeward  part  of  the  configuration  about  in  the  middle  of  the  body.  According  to 
Reference  27,  this  type  of  model  support  has  the  most  severe  adverse  impact  on  the  vehicle 
aerodynamics.  Aerodynamic  measurements  are  lower  for  normal  force  than  they  should  be  because 
the  strut  mount  adversely  impacts  the  pressure  in  the  leeward  plane  of  the  body,  causing  the  body 
to  lose  normal  force.  Figures  3 8  A,  38B,  and  38C  present  the  experimental  and  predicted  values  of 
Cn  and  Xcp  for  M,,  =  0.6,  0.9,  and  1.3,  respectively,  for  both  the  $  =  0  and  45  deg  roll  positions. 
Note  that  center  of  pressure  is  reasonably  well  predicted  at  all  Mach  numbers  and  at  both  0  =  0  and 
45  deg.  However,  normal  force  is  predicted  accurately  on  average  only  at  M  =  1.3.  The  average 
errorsonnormalforceforthethreeMachnumbersofO.6, 0.9,  and  1.3  are  19.3,  16.7,  and4percent. 
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FIGURE  36B.  SCHEMATIC  OF  REFERENCE  32  BODY-TAIL  CONFIGURATION 
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FIGURE  37.  NORMAL  FORCE  COEFFICIENT  AND  CENTER  OF 
PRESSURE  COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(CONFIGURATION  OF  FIGURE  36A;  M„  =  0.15) 
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FIGURE  38C.  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 
COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(CONFIGURATION  OF  FIGURE  36B;  M„  =  1.3) 
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respectively,  for  an  overall  average  for  these  three  Mach  numbers  of  13.3  percent.  The  average 
center  of  pressure  errors  at  all  three  Mach  numbers  is  a  quite  acceptable  2.6  percent  of  body  length 
and  is  fairly  consistent.  This  indicates  any  loss  in  normal  force  on  the  body  due  to  the  strut  at  high 
AOA  must  be  somewhat  evenly  distributed  along  the  body.  It  is  interesting  to  note  that  if  the 
subsonic,  high  AOA  («  >  55  deg)  normal  force  data  are  eliminated  from  the  averages,  the  overall 
normal  force  error  is  9  percent,  which  is  within  the  ±10  percent  tolerance  level. 

The  fourth  configuration  (Figure  39A)  is  a  wing-body  tested  at  both  NASA/LRC'*  and  New 
York  University.^®  A  sting  mount  was  used  and  maximum  AOA  was  50  deg.  Only  M„  =1.6  and 
2.7  data  were  available.  However,  oil  flow,  pressure  distribution  and  schlieren  data  were  also  taken. 
These  data  are  very  useful  in  showing  the  shock  structure  of  the  flow  field  and  in  particular  the 
interaction  of  the  bow  shock  wave  with  the  windward  plane  wings.  It  was  this  data  that  helped  the 
author  to  attempt  to  quantify  the  effects  of  "fin  choking"  at  the  O  =  45  deg  roll  position.  Figures 
39B  and  39C  present  the  comparison  of  predicted  normal  force  coefficients  and  centers  of  pressure 
at  M„  =1.6  and  2.7  respectively.  As  seen  in  the  figures,  excellent  agreement  with  the  data  is  shown 
for  both  the  <&  =  0  and  45  deg  roll  positions.  Average  errors  are  3.7  percent  and  0.6  percent  i  for 
Cn  and  X(,p  respectively.  It  is  also  interesting  to  note  that  there  is  very  little  difference  in  the 
aerodynamics  for  the  $  =  0  and  45  deg  plane  up  to  the  point  where  "fm  choking"  starts  to  occur,  at 
which  point  the  $  =  45  deg  Cn  data  are  slightly  lower  than  <&=  0,  and  the  Xc?  is  slightly  less  stable 
for  O  =  45  deg. 

The  configuration  of  Figure  39A  had  an  r/s  of  0.2,  which  meant  it  was  wing-dominated.  For 
wing-dominated  configurations,  the  difference  in  aerodynamics  between  the  <&  =  0  and  45  deg 
planes  is  much  smaller  than  for  configurations  that  have  large  values  of  r/s.  It  is  believed  that  is  in 
large  part  due  to  the  lower  minimum  value  of  Kb(w)  at  high  AOA  for  $  =  45  deg  versus  0  =  0  deg, 
as  illustrated  in  Figure  24. 


FIGURE  39A.  WING-BODY  CONFIGURATION  USED  IN  VALIDATION  PROCESS 
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FIGURE  39B.  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 
COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(M„=1.6) 


FIGURE  39C.  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 
COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(M„  =  2.7) 
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The  last  three  single-fm  configurations  are  shown  in  Figure  40A  and  consist  of  a  10-caliber 
body  with  a  3.5-caliber  tangent-ogive  nose  followed  by  a  6.5-caliber  afterbody.  Six  tail  planforms 
were  tested  in  Reference  33  and  three  of  those  are  illustrated  in  Figure  40 A.  The  three  tail 
planforms  allow  a  variation  in  aspect  ratio  of  0.31  to  1.24  with  r/s  varying  from  0.6  to  0.27 
respectively.  The  normal  force  coefficient  for  the  configurations  of  Figure  40A  are  shown  in 
Figure  40B  for  both  =  0  and  45  deg  roll  at  AOA  12  and  20  deg  and  as  a  function  of  Mach 
number.  As  noted  in  the  Figure  40B,  good  agreement  of  the  theory  and  experiment  is  shown  for 
both  AO  As,  at  both  roll  positions  and  for  the  Mach  number  range  of  1.5  to  4.6  where  data  were 
available.  It  is  also  worthy  of  note  that  the  O  =  0  and  <&  =  45  deg  roll  results  differ  primarily  for 
the  higher  r/s  and  lower  Mach  number.  This  is  due  primarily  to  the  lower  minimum  body  carryover 
lift  for  $  =  45  deg  compared  to  O  ==  0  deg  at  high  AOA.  Higher  AOA  data  than  20  deg  were  not 
available  for  this  configuration. 


3.2  WING  OR  CANARD-BODY-TAIL  CONFIGURATIONS 

Several  configurations  with  two  sets  of  lifting  surfaces,  with  and  without  control  deflections 
are  also  considered  for  the  validation  process.  The  first  of  these  is  an  old  version  of  the  standard 
missile  configuration  tested  at  M  =  0.1  by  Dunn,  et  al.^'*  at  the  Naval  Postgraduate  School. 
Reference  34  showed  only  normal  force  coefficient  data,  but  data  were  given  to  AOA  90  deg  at  both 
the  $  =  0  and  45  deg  roll  orientations.  Figure  41  shows  the  configuration  tested  in  the  wind  tunnel 
and  the  configuration  modeled  in  the  IAP95.  Experience  has  shown  that  when  a  wing  planform 
does  not  fit  the  specifications  of  the  IAP95  logic,  then  the  parameters  to  hold  constant  are  wing  area, 
span,  aspect  ratio,  centroid  of  planform  area  and  taper  ratio.  Other  parameters  (chord  length  and 
sweep  angles)  can  be  varied  slightly  to  meet  the  constraints.  Comparison  of  theoretical  and 
experimental  values  of  normal  force  coefficient  as  a  function  of  AOA  and  at  $  =  0  and  45  deg  roll 
are  shown  at  the  bottom  of  Figure  41 .  Once  again  the  comparisons  are  quite  good  up  to  AOA  of 
40  to  50  deg  where  the  theory  begins  to  be  lower  than  the  data  suggest.  However,  since  the  model 
was  sting-mounted,  it  is  suspected  that  sting  interference  effects  contribute  to  some  of  this  difference 
between  data  and  theory  at  the  high  AOA. 

It  is  also  interesting  to  note  the  fairly  large  difference  between  the  $  =  0  and  45  deg  normal 
force  coefficients  at  high  AOA.  This  difference  can  be  as  much  as  40  to  50  percent  and  is  attributed 
primarily  to  the  dorsal  body  carryover  lift,  which  gives  a  much  higher  value  at  $  =  0  versus  <&  =  45 
deg.  Figure  24,  for  large  r/s  values,  illustrates  this  qualitatively. 

A  second  wing-body-tail  case  is  shown  in  Figure  42.  Data  from  this  configuration  were 
available^^  only  for  0  =  0,  but  the  case  was  a  good  candidate  to  compare  the  aerodynamics  of  the 
AP95  to  the  IAP95,  due  to  its  having  fairly  large  wings  and  tails  that  were  fairly  close  together 
axially.  The  wing-tail  model  primarily  affects  pitching  moments,  as  the  actual  normal-force  losses 
on  the  tail  due  to  presence  of  the  forward  lifting  surface  can  be  fairly  small  in  percentage  except  at 
small  AOA.  As  a  result,  only  pitching  moments  are  shown.  These  are  given  in  the  lower  portion 
of  Figure  42  for  M»  of  1.41,  2.0,  2.5,  and  3.08.  These  Mach  numbers  bracketed  the  data  from 
Reference  35  that  were  taken  to  AOA  of  25  deg.  As  seen  in  Figure  42,  the  new  wing-tail  model 
gives  some  improvement  over  the  AP95  model,  with  the  improvements  being  the  most  noticeable 
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FIGURE  40A.  THREE  BODY-TAIL  CONFIGURATIONS  FROM  REFERENCE  33 
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FIGURE  40B.  NORMAL  FORCE  COEFFICIENT  VS.  MACH  NUMBER  FOR 
CONFIGURATION  OF  FIGURE  40A 
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CONFIGURATION  TESTED  IN  WIND  TUNNEL 
(FROM  REFERENCE  34  WHERE  DIMENSIONS  ARE  IN  INCHES) 


->]  1<-  0.26 


MODIFIED  CONFIGURATION  USED  IN  AEROPREDICTION  COMPUTATIONS 
PARAMETERS  FOR  BOTH  MODELS 

(AR)j=4.0  bt  =  3.76in.  Xj=0.16  (Ale)t=24°  AT=3.54in? 

(AR)q  =  0.12  bQ  =  1.32  In..  Xq  =  0.77  (Ale)d  =  ^D  = 


FIGURE  41.  NORMAL  FORCE  COEFFICIENT  COMPARISON  OF  THEORY 
AND  EXPERIMENT  (M„  =  0. 1) 
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as  AOA  increases  above  20  deg.  This  is  because  the  improved  model  discussed  in  section  2.2.8 
requires  the  wing-tail  interference  term  due  to  AOA  to  go  to  zero  much  faster  than  the  model  in  the 
AP95.  This  was  based  on  the  data  of  Reference  22,  which  is  seen  to  hold  for  the  configuration  of 
Reference  35  as  well. 

The  next  two  configurations  are  shown  in  Figure  43  with  the  wind  tunnel  data  given  in 
Reference  36.  Both  of  these  configurations  have  a  fairly  large  lifting  surface  of  aspect  ratio  close 
to  one.  The  configuration  in  Figure  43A  has  a  fairly  large  aspect  ratio  tail  for  control  whereas  the 
configuration  in  Figure  43B  has  a  moderate  aspect  ratio  canard  for  control.  Note  that  the  canard 
modelled  in  the  IAP95  is  slightly  different  than  that  tested  in  the  wind  tunnel  due  to  the  allowable 
variables  in  the  IAP95.  This  canard  has  the  same  area,  aspect  ratio,  span,  and  centroid  of  presented 
area  as  that  tested  in  the  wind  tunnel. 

The  comparison  of  predicted  lift  and  pitching  moment  coefficients  with  data^®  for  the 
configuration  of  Figure  43  A  is  given  in  Figure  44.  Figure  44 A  gives  the  Mach  1.6  and  2.0  results 
at  AOA  of  10  deg;  Figure  44B,  the  Mach  2.36  and  2.86  results  at  AOA  of  20  deg;  and  Figure  44C, 
the  Mach  3.95  and  4.63  results  also  at  AOA  20  deg.  The  results  are  shown  as  a  function  of  control 
deflection,  since  data  were  available  for  control  deflections  up  to  30  deg  for  all  conditions  tested, 
whereas  data  were  available  only  up  to  AOA  of  10  deg  at  lower  Mach  numbers  and  slightly  over 
20  deg  at  higher  AOA.  Hence,  both  Figure  43  configurations  serve  to  test  the  nonlinear  control 
deflection  prediction  model  with  some  AOA  nonlinearities  also  included  (due  to  an  upper  AOA  of 
only  20  deg).  As  seen  in  Figure  44,  at  all  Mach  numbers,  AO  As  and  6s,  both  Cl  and  Qi  are 
predicted  well  within  the  accuracy  goal  of  ±10  percent  on  normal  force  and  ±4  percent  of  body 
length  on  center  of  pressure.  While  the  pitching  moment  comparisons  in  the  figures  do  not  look 
quite  as  good  as  normal  force,  when  one  translates  this  to  center  of  pressure,  results  are  quite 
adequate.  The  worst-case  center  of  pressure  prediction  in  Figure  44  is  the  6  =  30  deg,  <&  =  0  deg, 
M  =  2.36,  and  AOA  =  20  deg  case  shown  in  Figure  44B.  For  this  case,  the  Xc?  error  is  only 
2.8  percent  of  body  length.  Since  the  body  length  of  both  configurations  in  Figure  43  is 
14.55  calibers,  the  worst  case  error  for  Xcp  is  0.4  caliber. 

The  other  point  worthy  of  note  in  Figure  44  is  the  comparison  of  the  $  =  45  deg 
aerodynamics  with  those  of  $  =  0  deg.  In  general,  the  lift  coefficient  is  lower  for  $  =  45  deg  than 
for  $  =  0  deg.  The  lower  the  Mach  number,  the  greater  the  difference  between  the  two  roll 
orientations.  This  is  due  to  a  combination  of:  lower  minimum  values  of  body  carryover  lift  for 
$  =  45  deg  than  for  $  =  0  (see  Figure  24);  lower  values  of  wing-tail  interference  at  high  Mach 
number  and  AOA;  and  large  control  surfaces.  In  essence,  the  physical  mechanisms  that  cause 
differences  in  aerodynamics  between  <&  =  0  and  45  deg  seem  to  dissipate  as  Mach  number  gets  into 
the  hypersonic  flow  regime,  where  the  Newtonian  Impact  Flow  assumption  appears  to  be 
reasonable.  Note  also  the  lower  stability  of  the  $  =  45  roll  aerodynamics  compared  to  the  =  0 
case  at  all  Mach  numbers  in  Figure  44. 

The  comparisons  of  theory  with  experiment  for  the  canard-control  configuration  of 
Figure  43B  are  given  in  Figure  45.  This  figure  is  similar  to  Figure  44  in  that  results  of  lift  and 
pitching  moment  coefficients  are  given  for  Mach  1.6,  2.0,  2.36,  2.06,  3.95,  and  4.63  for  the  0  =  0 
and  45  deg  roll  orientations,  then  compared  to  data.  Here  the  control  deflection  is  positive  (leading 
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1.17 


FIGURE  43 A.  WING-BODY-TAIL  CONFIGURATION  USED 
IN  VALIDATION  PROCESS^*  (DIMENSIONS  IN  INCHES) 


1.3 


FIGURE  43B.  CANARD- WING-BODY  CONFIGURATION  USED 
IN  VALIDATION  PROCESS^*  (DIMENSIONS  IN  INCHES) 
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FIGURE  45C.  LIFT  AND  PITCHING  MOMENT  COEFFICIENTS  FOR 
FIGURE  43B  CONFIGURATION  AT  AOA  =  20  DEG 
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edge  up)  rather  than  negative  as  in  Figure  44.  Also,  the  control  surfaces  have  a  much  smaller  area 
than  the  tail  surfaces  of  Figure  43  A.  These  slightly  smaller  control  surfaces  also  tend  to  decrease 
the  difference  between  the  <>  =  45  deg  and  $  =  0  deg  aerodynamics  due  to  lower  wing  tail 
interference  as  well  as  lower  forces  from  control  deflections.  Comparisons  of  theory  and  data  once 
again  are  quite  adequate  and  within  the  desired  accuracy  bounds.  The  worst  case  comparisons  are 
for  the  higher  Mach  number  lift  coefficient  predictions  and  large  6’s.  Notice  in  Figure  45C  that 
experimental  data  show  Cl  actually  decreasing  with  increasing  6  at  $  =  45  deg.  This  gives  errors 
in  predictions  of  up  to  10  percent  in  Cl-  It  is  believed  the  reason  for  this  phenomena  is  shock 
interference  from  the  forward  deflected  control  surfaces  onto  the  fairly  large  aft  wings.  This 
phenomena  is  not  accounted  for  in  the  IAP95. 

The  final  configuration  shown  is  the  SEASP ARROW  configuration  as  shown  in  Figure  46A 
with  wind  tunnel  results  taken  from  Reference  37.  This  configuration  has  fairly  large  wings  and 
tails,  with  the  wings  used  as  the  control.  Data  were  taken  at  Mach  numbers  of  1.5,  2.0,  2.35,  2.87, 
3.95  and  4.63.  Figures  46B,  46C,  and  46D  give  comparisons  of  the  IAP95  with  the  data  at  M  =  1.5, 
2.87,  and  4.63  respectively;  Results  are  shown  for  control  deflections  of  0  and  10  deg  in  Figure  46B 
and  for  0  and  20  deg  deflection  in  Figures  46C  and  46D.  Also,  all  results  for  both  the  $  =  0  and 
45  deg  roll  positions  are  shown.  Several  comments  are  in  order  with  respect  to  the  overall 
comparisons.  First  of  all,  the  LAP  model  achieves  its  goal  of  predicting  average  accuracy  of  C^,  Cjj 
of  ±10  percent  and  Xcp  of  ±4  percent  ^  on  this  configuration  at  both  the  $  =  0  and  45  deg  roll 
positions.  Secondly,  Cn  is  predicted  equally  well  at  $  =  0  and  45  deg.  However,  Cm  is  predicted 
better  at  $  =  45  deg  than  at  $  =  0  deg  due  to  the  center  of  pressure  shift  derived  in  section  2.2.6. 
No  such  shift  has  been  applied  at  O  =  0  deg.  Apparently  one  is  needed,  but  it  is  not  clear  what  the 
physical  justification  is.  As  seen  in  the  pitching  moment  predictions  for  the  $  =  0  deg  roll 
orientation,  the  IAP95  in  general  is  slightly  too  stable.  The  third  point  is  that  Ca  prediction  is 
slightly  better  for  the  O  =  0  deg  roll  orientation  than  the  <&  =  45  deg  position.  Apparently,  the  factor 
of  1.414  applied  to  the  fins  for  control  deflection  component  of  axial  force  is  quite  adequate  at  low 
angle  of  attack  at  all  Mach  numbers  as  seen  in  the  figures.  However,  at  high  angle  of  attack,  this 
factor  appears  to  be  too  high  at  the  lower  supersonic  Mach  numbers  and  too  low  at  the  higher 
supersonic  Mach  numbers.  It  is  suspected  that  at  the  lower  supersonic  Mach  numbers,  only  the  fins 
in  the  windward  plane  should  have  the  full  factor  of  yjl,  whereas  those  in  the  leeward  plane  should 
have  a  lower  factor.  For  high  supersonic  conditions,  it  is  suspected  the  bow  shock  and  internal 
shock  interactions  actually  add  to  the  factor  of  y'T.  An  empirical  model  for  Ca  to  account  for  this 
physics  would  improve  the  Ca  comparison,  but  time  does  not  permit  this  effort  at  present. 

The  final  point  to  be  made  is  concerning  the  high  Mach  number,  high  AOA  conditions  in 
Figure  46D.  Notice  that  the  pitching  moment  suddenly  loses  stability  above  AOA  of  30  deg  and 
$  =  0  deg  roll.  The  normal  force  also  decreases  somewhat  as  well.  It  is  believed  this  is  due  to 
internal  shock  waves  from  the  bow  shock  and  wing  shock  intersecting  the  tail,  causing  a  loss  of  tail 
lift  and  a  sudden  decrease  in  stability.  No  accounting  of  these  internal  shocks  from  a  forward 
mounted  fin  to  an  aft  mounted  fin  is  made  in  the  IAP95.  It  is  also  interesting  to  note  that  above 
AOA  70  deg,  while  the  results  are  not  shown,  the  IAP95  predictions  for  pitching  moment  agree 
quite  well  with  the  data.  Apparently  the  internal  shock  interactions  are  important  between  AOA  of 
about  30  to  60  deg. 
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4.0  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


New  technology  has  been  developed  to  allow  engineering  estimates  of  aerodynamics  of  most 
tactical  weapon  concepts  for  the  <&  =  45  deg  roll  position  (fins  oriented  in  "x"  or  cross-fm 
arrangement).  New  technology  developed  for  the  O  =  45  deg  roll  position  includes: 

a)  Nonlinear  wing-body  and  body-wing  interference  factor  methodology  due  to  AOA 
and  control  deflection. 

b)  A  new  semiempirical  wing-tail  interference  model  for  both  the  O  =  0  and  $  =  45 
deg  roll  positions. 

c)  An  approximate  method  to  estimate  wing  alone  shift  in  center  of  pressure  at  $  =  45 
deg  and  high  AOA. 

d)  An  approximate  method  to  estimate  "fm  choking"  at  high  M  and  AOA  based  on 
Newtonian  impact  theory 

This  new  technology,  along  with  the  $  =  0  methodology  of  Reference  1,  now  allows  trim 
aerodynamics  to  be  computed  in  the  two  roll-stabilized  planes  that  most  of  the  world’s  missiles  fly. 
The  data  bases  upon  which  the  methodology  was  based  was  limited  in  Mach  number  to  4.5  and 
AOA  of  40  deg.  However,  engineering  judgement  and  other  data  were  used  to  allow  calculations 
to  be  performed  to  AOA  90  deg,  Mach  numbers  up  to  20,  for  axisymmetric  solid  rocket  weapons 
with  up  to  two  sets  of  lifting  surfaces. 

Based  on  the  new  methodology  and  computations  to  date,  the  following  conclusions  are 

made. 


a)  The  two  primary  reasons  missile  normal  forces  are  lower  at  $  =  45  deg  than  at 
<&  =  0  deg  is  that  the  minimum  value  of  body  carryover  lift  at  high  AOA  is  lower  for 
<&  =  45  than  $  =  0.  The  larger  the  value  of  r/s,  the  more  difference  between  these 
minimum  values.  If  the  configuration  has  two  sets  of  lifting  surfaces,  the  wing-tail 
interference  is  higher  for  AOA  >  20  deg  for  the  $  =  45  deg  plane  than  for  the 
$  =  0  deg  plane,  also  contributing  to  a  lower  normal  force  for  O  =  45  deg  and  a  less 
stable  configuration. 

b)  At  high  AOA:  Kw(b)  approaches  1 .0  for  most  Mach  numbers;  Kb(W)  approaches  some 
minimum  value  that  is  a  function  of  r/s  and  Mach  number;  kw(B)  and  k^^^  are 
nonlinear  in  total  local  AOA  on  the  wing  and  are  functions  of  Mach  number  and 
total  AOA  on  the  wing. 

c)  The  actual  values  of  wing-tail  interference  are  much  smaller  than  those  predicted  by 
SBT  at  AOA  greater  than  about  20  deg,  particularly  for  Mach  numbers  greater  than 
about  2.0 
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d)  "Fin  choking"  appears  to  be  more  of  a  problem  on  larger  fm  (small  r/s) 
configurations  at  $  =  45  deg  than  on  smaller  fm  cases  (large  r/s). 

e)  For  Mach  numbers  less  than  one  and  AO  A  >30  deg,  it  is  not  clear  what  the  correct 
values  of  experimental  normal  force  are  for  a  given  configuration,  based  on  wind 
tunnel  results  available  in  the  literature. 

f)  Internal  shock  interaction  between  forward-mounted  and  aft-mounted  fins  becomes 
increasingly  important  as  both  AOA  and  Mach  number  increase.  The  current 
methodology  does  not  account  for  these  effects. 

g)  In  general,  ±10  percent  average  accuracy  has  been  maintained  for  both  normal  and 
axial  force  coefficients  and  4  percent  of  body  length  for  center  of  pressure  in  the 
$  =  45  deg  roll  position.  Exceptions  to  this  are  at  subsonic  Mach  number  and  high 
AOA  where  data  accuracy  is  in  question,  and  at  high  Mach  number  and  high  AOA 
for  configurations  that  have  two  sets  of  lifting  surfaces,  where  internal  shock 
interactions  may  be  a  problem. 

h)  The  current  overall  approach  of  using  linear  theory,  slender  body  theory,  or  second 
order  theory  for  low  angle  of  attack  aerodynamics,  and  estimating  the  nonlinear 
aerodynamic  terms  individually  and  directly  from  wind  tunnel  data  bases,  appears 
to  be  the  key  to  the  above  mentioned  average  accuracy  levels. 

Based  on  this  and  previous  research  (References  1,  5,  12,  and  16),  the  following 
recommendations  are  made  for  additional  work: 

a)  Additional  wind  tunnel  measurements  or  computational  fluid  dynamics  cases  need 
to  be  made  to  help  define  nonlinearities  of  the  interference  effects  as  a  function  of 
r/s. 

b)  A  method  is  needed  to  correct  wind  tunnel  data  at  subsonic  Mach  number  and  high 
AOA  for  sting  or  strut  mounting  effects. 

c)  An  engineering  method  is  needed  to  estimate  internal  shock  interaction  effects  for 
high  AOA  and  Mach  number. 

d)  Any  future  wind  tunnel  test  for  measuring  component  aerodynamics  should  be  done 
with  lifting  surfaces  large  enough  to  separate  out  body  and  wing  lift  accurately,  with 
wings  mounted  in  the  middle  of  the  body  and,  preferably,  with  simultaneous 
measurements  of  body  forces  in  conjunction  with  the  wing  and  wing  forces  in 
conjunction  with  the  body.  This  would  allow  more  accurate  determination  of  the 
interference  terms  directly,  without  subtraction  of  two  large  numbers  to  obtain  a 
small  term. 
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AOA 

APC 

AP81 

AP93 

AP95 

AR 

COP 

DOF 

IAP95 

LT 

MNT 

NASA/LRC 

NSWCDD 

NW 

SB,  SBT 

SW 

a 

3.4 

Ap 


6.0  SYMBOLS  AND  DEFINITIONS 

Angle  of  Attack 
Aeroprediction  Code 
Aeroprediction  1981 
Aeroprediction  1993 
Aeroprediction  1995 
Aspect  Ratio  =  bVA^ 

Center  of  pressure 

Degrees  of  freedom 

Improved  Aeroprediction  Code  1995 

Linear  Theory 

Modified  Newtonian  Theory 

National  Aeronautics  and  Space  Administration/Langley  Research  Center 
Naval  Surface  Warfare  Center,  Dahlgren  Division 
Number  of  wings 

Slender  Body,  Slender-body  Theory 

Single  wing 

Speed  of  Sound  (ft/sec) 

Constants  used  in  nonlinear  wing-alone  model 

Planform  area  of  the  body  in  the  crossflow  plane  (ft^) 
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'■REF 


'WETIED 


Aw,  Sw 
Av 

b 

Ca 


Ca  ,C^ 


Tj  r.j. 


-'M 


Cn 

c 


N„ 


N, 


BCV) 


r  c 

Nbp^)’ 


AC 


N, 


B(W) 


"N, 


'N, 


T(V) 


'N„ 


'N, 


'WCB) 


Reference  area  (maximum  cross-sectional  area  of  body,  if  a  body  is  present, 
or  planform  area  of  wing,  if  wing  alone)(ft^) 

Planform  area  of  wing  in  crossflow  plane  (ft^) 

Area  of  body  or  wing  which  flow  touches 

Wing  span  (not  including  body)(ft) 

Axial  force  coefficient 

Base,  skin-friction,  and  wave  components,  respectively,  of  axial  force 
coefficient 

Crossflow  drag  coefficient 

Laminar  and  turbulent  skin  friction  coefficients  respectively 
Roll  damping  moment  coefficient 

Pitching  moment  coefficient  (based  on  reference  area  and  body  diameter,  if 
body  present,  or  mean  aerodynamic  chord,  if  wing  alone) 

Normal  force  coefficient 

Normal  force  coefficient  of  body  alone 

Negative  afterbody  normal-force  coefficient  due  to  canard  or  wing-shed 
vortices 

Normal-force  coefficient  on  body  in  presence  of  wing  or  tail 
Additional  normal-force  coefficient  on  body  due  to  presence  of  the  wing 
Linear  component  of  normal-force  coefficient 
Nonlinear  component  of  normal-force  coefficient 

Negative  normal-force  coefficient  component  on  tail  due  to  wing  or  canard- 
shed  vortex 

Normal  force  coefficient  of  wing  alone 
Normal-force  coefficient  of  wing  in  presence  of  body 
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'No 


Normal-force  coefficient  derivative 


Cp 

C. 


Cr 

Ct 

cal 

d 

dref 

dK 


w(B) 


d« 


d« 


deg 


El,  E2 

F)  Cl,  C2,  C3 

f,h 
fw,  fr 


Kb(w))  Kb(x) 


kB(W),  ^8(1) 


[Kb(W)]  MIN 


Kw(B);  Kx(b) 


Pressure  coefficient  ^  P” 

Stagnation  pressure  coefficient 
Side  force  coefficient 
Root  chord  (ft) 

Tip  chord  (ft) 

Caliber(s)  (one  body  diameter) 

Body  diameter  (ft) 

Reference  body  diameter  (ft) 

Rate  at  which  Kw(B)  or  Kb(W)  decreases 
Degree(s) 

Constants  used  in  semiempirical  wing-tail  interference  model 

Dimensionless  empirical  factors  used  in  nonlinear  models  of  and 
to  approximate  effects  due  to  high  AO  A  or  control  deflection. 

Lateral  and  vertical  position  of  wing  vortex 

Lateral  location  of  wing  or  tail  vortex  (measured  in  feet  from  body  center  line) 
Tail  interference  factor 

Ratio  of  additional  body  normal-force  coefficient  in  presence  of  wing,  or  tail 
to  wing,  or  tail  alone  normal-force  coefficient  at  6  =  0  deg 

Ratio  of  additional  body  normal-force  coefficient  due  to  presence  of  wing  or 
tail  at  a  control  deflection  to  that  of  wing  or  tail  alone  at «  =  0  deg 

Minimum  value  of  Kb(w)  as  percent  of  slender-body  theory  value 

Lowest  value  of  [Kb(w)]min  for  O  =  45  deg 

Ratio  of  normal-force  coefficient  of  wing  or  tail  in  presence  of  body  to  that 
of  wing  or  tail  alone  at  6  =  0  deg 
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kw(B);  kx(B) 

[AK-vv(B)]„=0 

and 

[AKb(w)]„=o 

C 

«N 

m 

M 

Me 

Ml 

P 

Pjj  Pw 

Po,  pL 

r 

% 

fwj  Tx 

Re 

R^CRIT 

Rcd 

R^eff 

s 

t/Cr 


Ratio  of  wing  or  tail  normal-force  coefficient  in  presence  of  body  due  to  a 
control  deflection  to  that  of  wing  or  tail  alone  at «  =  0  deg 

Amount  that  the  experimental  values  of  Kw(b)  and  Kb(W)  exceed  slender  body 
theory  at «  =  0“ 

Length  (ft) 

Nose  Length  (can  be  in  calibers  or  feet) 

Cot  Ale 

Mach  number  =  V/a 

Normal  Mach  number  to  body  axis  =  M  sin  « 

Local  Mach  number 
Pressure  (Ib/ff ) 

Loading  factors  in  leeward  and  windward  planes  respectively 

Local  and  stagnation  pressure  respectively  (Ib/ff ) 

Radius  of  body  (ft) 

Radius  of  nose  tip  (ft) 

Radius  of  body  at  wing  or  tail  locations 
pVC 

Reynolds  number  =  - — 

Critical  Reynolds  number  where  flow  transitions  from  separating  on  the 
forward  part  of  a  circular  cylinder  to  the  rear  part 

Reynolds  number  based  on  diameter  of  body 

Reynolds  number  defined  to  better  correlate  values  of  versus  Reynolds 
number 

Wing  or  tail  semispan  plus  the  body  radius  in  wing-body  lift  methodology 
Tail  thickness  to  its  root  chord 
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t/d 

TW 

U,  V,  w 

V 

X 

XcP 


(Xcp)l)  (Xcp)nl 

x,y,2 


Y‘,  Z' 

YcP 


“c 

“*0 


“f 

“n 


®'W)  “'x 


“b  “2 


P 

6 


Tail  thickness  to  body  diameter 
Bevel  length  on  leading  edge  of  wind  tunnel  model 
Perturbation  velocities  in  x,  y,  z  directions 
Velocity  (ft/sec) 

Distance  along  the  axis  of  symmetry  measured  positive  aft  of  nose  tip  (ft 
or  cal) 

Center  of  pressure  (in  feet  or  calibers  from  some  reference  point  that  can  be 
specified)  in  x  direction 

Center  of  pressure  of  linear  and  nonlinear  terms  of  normal  force. 

Axis  system  fixed  with  x  along  centerline  of  body 
Body  fixed  axis  system  with  X  ■  along  centerline  of  body 
Center  of  pressure  of  wing  lift  in  y  direction 
Angle  of  attack  (deg) 

Angle  of  attack  where  wing-body  interference  factor  starts  decreasing  (deg) 

Angle  of  attack  where  the  wing-body  interference  factor  reaches  a  minimum 
(deg) 

Value  of  «  (as  percent  of  «n)  where  Cj^  reaches  a  maximum  magnitude 
Angle  of  attack  where  Kw(B)  reaches  a  constant  value 
Value  of  «  where  C„  goes  to  zero 

^T(V) 

Value  of  «  where  goes  to  zero  for  At^/Aref  =  5.5 

^T(V) 

Local  angle  of  attack  of  wing  or  tail  (“^^  +  8  or  °=.j,  +6,  respectively,  in 
degrees) 

Angles  of  attack  used  in  nonlinear  model  for  Kb(w) 

-  1 

Control  deflection  (deg),  positive  leading  edge  up 
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6eq 


67 


n 


Y 

r 

k 

■^LE 

Mo 

$ 

<l> 

p 

0 


00 


q. 


Angle  between  velocity  vector  and  tangent  to  body  surface 

Deflection  of  wing  or  tail  surfaces  (deg),  positive  leading  edge  up 

Parameter  used  in  viscous  crossflow  theory  for  nonlinear  body  normal  force 
(in  this  context,  it  is  the  normal  force  of  a  circular  cylinder  of  given  length-to- 
diameter  ratio  to  that  of  a  cylinder  of  infinite  length) 

Specific  heat  ratio 

Vortex  circulation,  positive  counter  clockwise  facing  upstream  (ff/sec) 
Taper  ratio  of  a  lifting  surface  =  cjc^ 

Leading  edge  sweep  angle  of  wing  or  tail  (deg) 

Mach  angle  =  sin'*  (1/M) 

Viscosity  of  air  (slug/ft-sec) 

Circumferential  position  around  body  where  $  =  0  is  leeward  plane  (deg) 
Velocity  potential 
Density  of  air  (slug/ff ) 

Local  surface  slope  of  body  with  respect  to  body  axis 
Free-stream  conditions 

Dynamic  pressure  of  freestream  =  I/2  p.  V„^  (Ib/ft^) 
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APPENDIX  A 

SLENDER  BODY  AND  LINEAR  THEORY  RESULTS 
FOR  ROLL-DEPENDENT  AERODYNAMICS 


A-l/A-2 
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For  the  slender  body  and  linear  theory  analysis,  References  A-1  and  A-2  will  be  used.  The 
wing-alone  and  wing-body  theories  will  be  examined  separately. 


WING-ALONE  AERODYNAMICS 

(^N  ^6  of  Equation  (2)  in  the  main  text  are  all  normal  force 

coefficients  of  the  wing  or  tail  in  a  uniform  freestream  with  no  interference  effects.  Here  the 
subscripts  «  and  6  refer  to  normal  force  coefficient  due  to  AOA  or  control  deflection.  These  terms 
all  have  a  linear  and  a  nonlinear  component.  Disregarding  the  subscripts  to  distinguish  wing  from 
tail  or  angle  of  attack  from  control  deflection,  the  normal  force  of  the  wing  can  be  generically 
defined  as  consisting  of  a  linear  and  a  nonlinear  term.  That  is. 


Cn  =  C. 


+  c. 


(A-1) 


The  linear  term  of  Equation  (A-1)  is  normally  computed  from  integrating  the  pressure  differential 
between  the  lower  and  upper  part  of  the  wing  surface  over  the  entire  wing.  It  is  instructive  to  review 
some  of  the  fundamentals  of  thin  wing  theory  to  help  understand  how  roll  orientation  influences  (or 
does  not  influence)  the  linear  term  of  Equation  (A-1). 

The  mathematical  model  to  determine  Cj^  is  based  on  small  perturbation  theory  for  both 
subsonic  and  supersonic  flow.  In  subsonic  flow,  solutions  to  the  perturbation  equation  are  referred 
to  as  lifting  surface  theory,  whereas  in  supersonic  flow  they  are  referred  to  as  three-dimensional  thin 
wing  theory.  The  small  perturbation  equation  of  motion  is 


(l  -  Mf)  <1>„  +  <1>„  +  <1>„  =  0 


(A-2) 


The  boundary  conditions  require  the  flow  to  be  tangent  to  the  surface, 

=  <1^2  =  [  — ]  +  <="  +  dg/dX  (A-3A) 


and  the  perturbation  potential  and  velocities  to  vanish  at  a  large  distance  from  the  wing. 

(1)  =  (j)^  =  cf)^  =  (})^  =  0  at  x,y,z  =  (A-3B) 


A-3 
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The  flow  tangency  boundary  condition,  Equation  (A-3a),  can  be  separated  into  terms  due  to 

thickness,  chamber,  and  angle  of  attack.  For  most  missiles,  the  airfoil  is  symmetric  about  the  x-y 

plane  so  the  chamber  effect  is  zero.  The  normal  approach  to  solution  of  Equation  (A-2)  is  to 

solve  the  thickness  —  |  and  lift  problems  («)  independently,  since  Equation  (A-2)  is  linear.  Then, 
dx  / 

if  desired,  the  solutions  can  be  linearly  added  together. 

The  linearized  pressure  coefficient  for  a  flat  surface  from  Equation  (A-2)  is 

Cp  =  -2  (|)^  =  -2u(x,y,z)  (A-4) 

For  the  axial  force  problem,  the  boundary  condition  is  symmetric  about  the  plane  z  =  0.  For  the  lift 
problem,  the  normal  force  is  the  difference  between  the  pressure  on  the  upper  and  lower  surfaces 
of  the  wing.  For  this  case 


ACp  =  4«  f  (M„,  x,y,A)  (A-5) 

The  function  f  (M„,  x,  y.  A)  is  used  here  to  indicate  that  a  complicated  expression  exists  to  define 
the  differential  pressure  from  upper  to  lower  wing  surfaces  that  is  a  function  of  M,,  and  the  wing 
geometric  parameters.  This  function  is  different  for  subsonic  and  supersonic  leading  edge  wings. 

The  important  point  is  that  for  both  the  axial  and  normal  force  of  a  wing  in  isolation,  there 
is  no  roll  orientation  or  sideslip  dependence  in  the  pressure  coefficient.  As  a  result,  we  can  analyze 
wings  in  isolation  directly  without  concern  for  crosscoupling  effects. 

Examining  a  wing  (see  Figure  A-1  A)  at  some  roll  orientation,  it  is  found  that  in  the  direction 
normal  to  the  planform,  the  velocity  is  V„  sin  «  cos  and  angle  of  attack  is  sin  «  cos  or  for  the 
small  angle  approximation  of  sin  «  «  the  velocity  is  V„ cos  Also,  to  rotate  the  normal  force 
vector  (normal  to  the  lifting  surface  area)  to  the  vertical  plane,  another  cos  $  is  required. 

The  normal  force  of  the  wing  alone  in  the  x,y  plane  is  therefore 


(A-6) 


If  there  is  a  wing  perpendicular  to  the  planform  in  Figure  A-1  A,  so  that  a  cruciform  arrangement 
is  present,  then  for  this  wing  (see  Figure  A- IB) 
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FIGURE  A-IA.  LT  AND  SET  IMPLICATIONS  FOR  ROLL -DEPENDENT 
AERODYNAMICS,  SINGLE  WING  OF  SPAN  B 
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FIGURE  A- IB.  LT  AND  SBT  IMPLICATIONS  OF  ROLL-DEPENDENT 
AERODYNAMICS,  CRUICIFORM  WINGS  OF  SPAN  B 
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For  a  cruciform  missile,  the  Equations  (A-6)  and  (A-7)  are  added,  yielding 


or 


(A-8) 


This  result  says  that  for  small  angles  of  attack  for  a  cruciform  missile,  the  lift  of  the  wing  is 
independent  of  roll  angle.  This  is  also  true  of  the  axial  force  as  well.  This  means  that  we  can 
analyze  the  cruciform  missile  fins  alone  at  roll  orientation  of  zero  and  apply  these  results  to  any  roll 
orientation. 

It  is  also  worthy  to  note  in  Figures  (A-IA)  and  (A-IB)  the  implications  of  LT  and  SBT  on 
the  other  aerodynamics  for  wings  alone.  Of  particular  note  is  the  fact  that  for  a  single  planar  wing, 
a  side  force  is  generated  for  O  =  0  (Figure  A-IA)  whereas  for  a  planar  set  of  cruciform  wings,  no 
side  force  is  generated.  Also,  for  more  than  one  wing,  the  axial  force  is  simply  the  product  of  the 
value  for  a  single  wing  times  the  number  of  wings  present.  Slender  body  theory  also  gives  the  result 
that  the  roll  damping  coefficient  of  a  cruciform  arrangement  is  only  1 .62  times  the  single  fin  value 
versus  twice  the  value."^"' 

Figure  A- 1C  presents  the  results  for  three  sets  of  wings,  where  wing-to-wing  interference 
has  been  ignored.  Note  once  again  the  fact  that  the  aerodynamics  are  independent  of  roll.  Also  the 
aerodynamics  are  all  a  function  of  the  number  of  wings,  NW. 

The  LT  and  SBT  results  of  Figure  A-1  are  true  for  flowfields  where  perturbations  in  the  flow 
due  to  the  presence  of  the  wing  are  small  and  where  any  shock  waves  are  weak.  As  a  result  of  these 
assumptions,  it  is  suspected  that  the  roll  independence,  roll  damping  and  other  force  results  of 
Figures  A- IB  and  A- 1C  may  not  hold  at  higher  angle  of  attack  where  the  aerodynamics  are  strongly 
nonlinear  and  fin-to-fm  interference  effects  are  higher.  On  the  other  hand,  this  assumption  will  be 
required  in  the  development  of  a  semiempirical  model.  The  reason  for  this  goes  back  to  the  way 
the  data  bases  of  References  A-4  to  A-6  were  obtained.  These  data  bases  obtained  aerodynamics 
on  a  single  wing  alone  at  O  =  0  and  wings  in  conjunction  with  the  body  at  various  roll  orientations. 
The  wing  in  conjunction  with  the  body  therefore  has  carryover  effects  from  the  body,  nonlinear  and 
linear  angle  of  attack  effects,  body  vortices,  fin-to-fin  interference,  and  any  roll  dependence  effects. 
Since  the  data  were  taken  in  a  total  sense,  it  would  be  difficult,  if  not  impossible,  to  accurately 
separate  all  the  independent  effects.  As  a  result,  a  quite  logical  way  to  use  the  Reference  6-9  data 
bases  is  through  using  the  wing-alone  data  at  O  =  0  for  any  roll  position,  combined  with  the  fin  in 
conjunction  with  the  body  data,  which  has  included  all  the  effects  mentioned  above.  The  net  result 
of  this  approach  will  be  to  include  all  the  nonlinear  effects  mentioned  into  the  interference  factors 
of  Equation  (2)  or  Equation  (3)  of  the  main  text. 
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0  =  0° 


© 


NW  =  3 


(Xcp)  ^0 

(^cp)  ~(^cp) 

V  '  '^/<D  =  o 

(Cy]  =0 

(Cy)  =0 

V  ^'(P=0 

FOR  MORE  THAN  3  SETS  OF  WINGS,  ALL  ABOVE  FORMULAS  ARE 

STILL  TRUE  EXCEPT  Ci  .  FOR  NW= 

'P 

(^'p)nw=co  "  '‘(^'p)sw 

FIGURE  A-IC  LT  AND  SET  IMPLICATIONS  FOR  ROLL-DEPENDENT 
AERODYNAMICS,  THREE  SETS  OF  WINGS  OF  SPAN  B 
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WING-BODY  AERODYNAMICS 

The  implications  of  slender  body  theory  on  the  wing-body  combination  will  now  be 
examined.  Particular  emphasis  will  be  placed  on  the  crossflow  plane.  The  small  perturbation 
equation  used  here  is  the  same  as  Equation  (A-2);  however,  the  boundary  condition  (A-3A)  is 
different.  The  reason  for  the  different  boundary  condition,  as  shown  in  Reference  A-2,  is  that 
in  the  perturbation  expansion  of  the  inner  and  outer  flow,  all  terms  to  the  first  order  in  the 
perturbation  parameter  e  vanish.  Thus  the  perturbation  velocities  in  the  outer  flow  are  of  the  order 
8^  versus  e  for  the  2-D  and  3-D  wing  cases.  Physically,  one  can  interpret  this  as  the  perturbation 
velocities  being  smaller  for  a  body  of  revolution  in  the  normal  force  direction  than  for  a  wing 
planform  because  the  body  allows  the  flow  to  wrap  around  it  in  a  smoother  way  than  a  wing.  This 
process  allows  the  perturbations  in  the  flow  at  a  large  distance  from  the  body  to  be  smaller  for  a 
body  than  for  a  wing. 

As  a  result  of  the  perturbations  being  of  the  order  versus  e,  the  pressure  coefficient 
equation  is  no  longer  linear  like  the  wing  alone.  The  approximate  pressure  coefficient  equation  for 
a  wing-body  combination  is  therefore 


Cp  =  -2(1)^  - 


=  -2u  -  (A-9) 

Recall  the  boundary  condition.  Equation  (A-3  A)  along  with  the  pressure  coefficient  Equation  (A-5), 
allowed  us  to  uncouple  the  thickness  boundary  condition  from  the  boundary  condition  due  to  angle 
of  attack.  Hence,  one  could  compute  axial  force  coefficients  with  Equation  (A-4)  with  only  the  first 
part  of  Equation  (A-3  A)  used  and  normal  force  coefficients  with  Equation  (A-5)  with  the  second 
part  of  Equation  (A-3A)  used.  This  can  no  longer  be  done  with  Equation  (A-9),  since  terms 
involving  the  axial  force  and  normal  force  boundary  condition  appear  to  the  second  power,  coupling 
the  axial  and  crossflow  solutions  together. 

A  fortunate  feature  of  triform  and  cruciform  missile  configurations,  is  that  when  the  force 
is  integrated  over  the  entire  configuration  in  the  lift  direction,  the  coupling  terms  cancel  out  and  the 
normal  force  is  thus  independent  of  roll  position^'*.  This  means  that  if  one  is  interested  in 
calculating  the  normal  force  of  a  wing-body  configuration  with  a  single  wing  (Figure  A-2  A),  the 
normal  force,  side  force  and  rolling  moment  are  all  roll-dependent.  However,  if  one  is  interested 
in  calculating  the  aerodynamics  of  a  cruciform  wing-body  combination  with  SBT  or  LT 
(Figure  A-2B),  the  aerodynamics  are  independent  of  roll. 

It  should  be  re-emphasized  that  as  long  as  trim  aerodynamic  models  are  of  interest  for 
missiles  flying  at  $  =  0  or  45  deg  roll  orientation,  and  symmetrical  deflections  of  the  control 
surfaces  are  considered,  the  approach  of  neglecting  the  out  of  plane  coupling  aerodynamics  is 
satisfactory.  However,  if  one  is  interested  in  attempting  to  use  an  engineering  code  to  generate 
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•  Cm  =  [Kw(B)  +  Kb(W)]  g  •  Cn  =  IKw{B)  +  Kb(W)]  cos^d)  +  Cn^ 

•  Ca  =  Ca  +Ca  •  Ca  =  Ca  +Ca 

«  «W  ''b  W  b 

•  =  Cy  =  0  •  Cy  =  [Kw(B)  +  Kb(W)1  cosO 

•  LEFT  AND  RIGHT  PANEL  LOADING  EQUAL  *  =  '  I^WCB)  Vcp  cos‘1’ 

w 


LOADING  ON  RIGHT  PANEL  GREATER 
THAN  THAT  ON  LEFT  PANEL  BY  AN  EQUAL 
AMOUNT  DUE  TO  OUT-OF-PLANE 
AERODYNAMIC  COUPLING 


FIGURE  A-2A.  SET  IMPLICATIONS  FOR  ROLL-DEPENDENT  AERODYNAMICS, 
WING-BODY  COMBINATION  (SINGLE  WING  OF  SPAN  B) 
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<I>=0.  8=0 


i  I 

Vc„sin  a 

•  Cn  =  [Kw(B)  +  Kb(W)]  Cn..,  +  Cn 

w  B 

• 

W  B 

•  =  Cy  =  0 

•  LEFT  AND  RIGHT  PANEL  LOADING  EQUAL 


8=0 

z 


VooSin  a 

•  Cn  =  [Kw(b)  +  Kb(w)1  Cn... + Cn 

W  B 

•  Ca  =  2Ca  +  Ca 

W  B 

•  Cy  =  =  0 

•  LOADING  ON  WINDWARD  PANELS 
GREATER  BY  AND  EQUAL  AMOUNT  THAN 
THAT  ON  LEEWARD  PANELS.  TERMS 
CANCEL  SOCy  =  C^  =  0 

• 

(CN)(I);>tO  =  (Cn)0=O 


FIGURE  A-2B.  SET  IMPLICATIONS  FOR  ROLL-DEPENDENT  AERODYNAMICS, 
WING-BODY  COMBINATION  (CRUCIFORM  WINGS  OF  SPAN  B) 
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aerodynamics  for  a  six-degree-of-freedom  trajectory  simulation  model,  the  out  of  plane 
aerodynamics  cannot  be  neglected.  An  approach  similar  to  that  in  Reference  A-7  will  be  required 
or  a  more  accurate  aerodynamic  approach  where  the  small  perturbation  assumptions  have  not  been 
made  will  most  likely  be  required. 


WING-BODY-TAIL  AERODYNAMICS  AT  $  -  0  DEG 

The  major  new  complexity  that  results  from  adding  two  sets  of  lifting  surfaces  to  a  missile 
configuration  (versus  one)  is  the  vortices  shed  from  the  forward  surfaces  impacting  the  rear  surfaces. 
For  this  discussion,  wings  and  canards  will  be  used  interchangeably  as  the  forward  surfaces,  and 
tails  will  be  the  aft  lifting  surfaces.  If  the  missile  is  at  supersonic  Mach  numbers  and  at  angle  of 
attack  greater  than  about  25  to  30  deg,  shock  wave  interference  between  the  forward  and  aft  lifting 
surfaces  or  between  the  bow  shock  wave  and  the  lifting  surfaces  can  also  occur.  This  problem  will 
not  be  discussed  in  the  linear  theory  and  slender  body  theory  discussion  but  is  discussed  in  the 
nonlinear  methods  discussion  of  the  report. 

If  the  missile  is  in  the  $  =  0  deg  roll  orientation,  only  one  wing  planform  will  release  a 
vortex  at «  >  0  deg  for  fairly  thin  wings  of  zero  chamber.  However,  at  ^  =  45  deg  for  a  cruciform 
wing  arrangement,  both  wings  release  vortices  and  these  vortices  follow  different  paths  depending 
on  whether  they  are  released  in  the  windward  or  leeward  plane.  It  is  instructive  to  examine  the  SBT 
results  for  both  the  $  =  0  and  45  deg  plane  in  terms  of  the  loss  of  lift  on  the  tail,  or  term  of 
Equation  (1).  In  reviewing  the  SBT  results,  References  A-1  and  A-8  will  be  the  primly  sources 
of  material. 

Figure  A-3,  taken  from  Reference  A-8,  is  a  schematic  of  the  $  =  0  deg  vortex  model  for  a 
planar  wing-body  combination.  The  wing-shed  vortices  are  shown  above  the  tail-body  cross  section. 
It  is  assumed  here  that  one  is  looking  upstream  of  a  wing-body-tail  configuration  at  a  positive  angle 
of  attack.  The  strength  of  the  wing-shed  vortices,  F,  is  assumed  to  be  positive  in  a  counterclockwise 
direction  looking  upstream. 

The  term  C„  of  Equation  (2),  is  defined  by'^’^ 

Mx(V) 


(A-IO) 


Some  of  the  parameters  of  Equation  (A- 10)  are  shown  in  Figure  (A-3).  The  parameter /n^,  the  point 
on  the  wing  where  the  single  vortex  separates,  is  approximated  by  the  centroid  of  the  approximate 
elliptic  circulation  distribution.  That  is 

fw  =  +  (^/4)  (bw/2)  (A-11) 


A-12 


NSWCDD/TR-95/160 


-r 


r 


(a) 


(b) 


(a)  WING  VORTICES  IN  CROSSFLOW  PLANE  OF  TAIL. 

(b)  TAIL  PLANFORM  DIMENSIONS. 


FIGURE  A-3.  MODEL  AND  DIMENSIONS  FOR  DETERMINATION  OF  TAIL  INTERFERENCE 
FACTOR  BY  STRIP  THEORY  FOR  4)  =  0  DEG  (FROM  REFERENCE  A-8) 
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Utilizing  Figure  A-4,  the  height  of  the  wing-shed  vortex  above  the  tail  center  of  area  is 


h  = 


°w 

2  J 


t“(4E>W  "  C,  [1 


j(l  -  x^han  -  -  gin  (6^)  ^(C^^/2 


tan  (Ale)w  + 


(C.) 


w 


b^/2 


'l] 

fO 

[4] 

.  2  , 

(A-12) 


Equation  (A-12)  assumes  the  vortex  leaves  the  wing  at  the  trailing  edge  outboard  location  given  by 
Equation  (A-1 1).  The  vortex  then  follows  the  velocity  vector  straight  back,  and  impacts  at  the  tail 
center  of  pressure  (here  approximated  along  the  root  mid  chord).  To  satisfy  the  boundary  condition 
on  the  body  surface  of  zero  normal  velocity  requires  that  the  vortex  of  equal  and  opposite  strength 
to  the  wing  vortex  be  placed  on  the  radius  vector  along  the  body  centroid  to  the  wing  vortex  and  at 
a  distance 


from  the  body  centroid.  Using  proportionality  of  triangles,  the  coordinates  of  the  image  vortex  are 


f  =  f 


f;  + 


(A-13) 


h,  =  h 


f;  + 


(A-14) 


A- 14 


VORTEX  PATH 


FIGURE  A-4.  NOMENCLATURE  FOR  DETERMINING  HEIGHT  OF  WING 
VORTEX  ABOVE  TAIL  MID  CHORD  FOR  «» =  0  DEG 
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It  should  be  noted  that  the  image  vortices  are  required  to  account  for  the  body  part  of  the 
configuration  in  the  vicinity  of  the  wing  and  are  equal  and  opposite  in  strength  to  the  wing-shed 
vortices. 

The  only  term  not  defined  in  Equation  (A- 10)  is  the  tail  interference  factor  if.  The  tail 
interference  factor  is  a  nondimensional  measure  of  the  interference  of  a  vortex  on  a  lifting  surface. 
In  mathematical  terms,  this  definition  is 


L  ) 
^T(V) 

(s,p  Tij.)  “ 

i  r 

(A-15) 


In  Equation  (A-15),  should  be  interpreted  as  the  lift  or  normal  force  on  the  tail-body  region  due 
to  a  vortex  shed  upstream;  should  be  interpreted  as  the  lift  on  the  tail  in  isolation.  Since  the 
velocity  of  the  vortices  is  counterclockwise  in  the  tail  region,  this  means  a  downward  velocity,  and 
hence,  the  lift  on  the  tail  due  to  the  wing-shed  vortex  is  negative.  Hence,  if  is  negative  and 
Equation  (A- 10)  is  negative. 

The  general  procedure  to  compute  L.j.^  is  through  the  use  of  strip  theory  since  this 
simpUfies  the  integration  of  the  lift  over  the  wing  surface  and  allows  for  a  closed  form  solution.  The 
lift  due  to  the  right  external  vortex  on  the  right  external  tail  panel  (looking  upstream  in  Figure  A-3) 
is 


L 


1 


4q„  r  (c,)t 


(A-16) 


The  function  L  of  Equation  (A-16)  is  defined  by 


(  r  f  1.  ) 

.  *x  n 

[(Sx  r.j,  A.J,)  (1  A.J.) 

+  (f^  -  B^)^' 

^  ®T  ®T  ®Ty 
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The  lift  on  the  right  panel  due  to  the  left  vortex  is  then 

4q„r  (c,) 


La  = 


/ 


2II  Wm:  -I  V. 


^  L 

/ 

A  —  - 

h  ) 
“2 

■)v„ 

V  ®T 

9 

®T 

®T,, 

(A-18) 


The  lifts  of  the  right  and  left  image  vortices  on  the  right  wing  panel  are  then  given  by 


L3  =  - 


4q.r  (c,)t  ^ 

fA 

211  (/m.=  -1)  V. 

9  9 

®T  ®T  ®T,i 

(A-19) 
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Since  the  left  tail  panel  has  the  same  lift  as  the  right  (because  of  symmetry  at  $  =  0  deg)  the  total 
lift  on  the  tail  due  to  the  vortices  shed  in  the  wing-body  region  is 


(A-21) 


Equation  (A-21)  abbreviates  the  parameters  in  the  L  function  Equation  (A- 17)  given  by  Equations 
(A- 16),  (A-18),  (A-19),  (A-20)  for  simplicity. 

The  lift  on  the  tail,  again  using  strip  theory  is 


4°=  q„o  (St  -  r  )  c  (1  +  A  ) 

T  —  *  i  ij  1 

-  1 


The  nondimensional  interference  factor,  ij,  then,  becomes 


It 


2 
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(A-23) 


A-17 


NSWCDD/TR-95/160 


It  should  be  noted  that  while  strip  theory  has  been  used  for  the  integration  of  the  lift  on  the  tail  in 
isolation  and  in  the  presence  of  the  vortex,  these  errors  tend  to  cancel  since  they  are 
nondimensionalized  by  one  another.  The  same  process  for  obtaining  the  interference  factors  Kw(b), 
Kb(W),  kw(B)  and  kB(W)  has  been  used  with  demonstrated  success. 

WING-BODY-TAIL  AERODYNAMICS  AT  O  =  45  DEG 

A  cruciform,  in-line,  wing-body-tail  arrangement  will  be  assumed  for  this  discussion.  This 
assumption  is  actually  not  necessary  for  the  slender  body  theory,  but  is  convenient  for  discussion 
purposes.  Figure  A-5  gives  the  geometrical  arrangement  for  vortex  placement  of  the  eight  vortices 
in  the  wing-body  plane,  which  must  be  accounted  for  in  the  tail  plane.  It  is  interesting  to  note  that 
even  though  slender  body  theory  gives  a  larger  lift  on  the  windward  plane  fins  than  on  the  leeward 
plane  fins  (see  wing-body  aerodynamics  discussion),  the  normal  approach  to  the  computation  of 
vortex  lift  on  the  tail  for  small  «  is  to  assume  (referring  to  Figure  A-5) 

r  =  -  r 


r,  =  -  r.  (A-24) 

This  being  the  case,  the  only  difference  between  the  $  =  0  deg  case  of  the  previous  section  and  the 
$  =  45  deg  case  is  the  difference  in  trajectories  due  to  the  roll  position.  The  physical  reason  that 
this  occurs  here  and  not  with  the  loading  computations  is  that  the  vortex  strength  is  a  function  of  the 
velocity  potential  difference  from  the  upper  to  the  lower  part  of  the  wing.  Those  potential 
differences  can  be  added  linearly  together.  On  the  other  hand,  the  section  2. 1 .2  results  showed  the 
velocity  to  the  second  power,  thereby  coupling  the  axial  and  crossflow  solutions  together.  It  is 
suspected  that  the  traditional  mathematical  approach  of  linearly  adding  the  velocity  potential 
differences  together  to  obtain  circulation,  is  one  reason  the  work  of  Reference  A-3  required  an 
empirical  correction  to  slender  body  theory,  particularly  as  angle  of  attack  increases.  This  is 
addressed  further  in  the  nonlinear  methods  discussion  of  the  report. 

For  the  $  =  45  deg  computation,  the  vortices  will  again  be  assumed  to  stream  rearward  along 
the  velocity  vector.  In  actuality,  the  windward  plane  vortices  actually  rise  above  this  line  because 
of  the  spin  direction;  the  leeward  plane  vortices  move  below  the  velocity  vector,  since  their  spin  is 
downward.  Reference  A-8  showed  that  the  more  sophisticated  method  of  vortex  tracking  was  not 
superior  to  the  simplified  approach,  however. 

Figure  A-5B  shows  the  vortex  pattern  of  the  eight  vortices  in  the  tail  plane  area  based  on  the 
simplified  assumption  that  their  paths  follow  the  velocity  vector.  It  is  of  interest  to  calculate  the 
effect  on  the  tail  surfaces  of  these  eight  vortices.  Since  the  vortex  pattern  is  symmetric  about  the 
x-z  plane,  only  the  tail  surfaces  in  the  right  half  plane  need  to  be  considered  and  these  results  then 
doubled  to  account  for  the  fins  in  the  left  half  plane  of  Figure  A-5B.  In  the  analysis,  since  the 
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FIGURE  A-5A.  VORTEX  MODEL  OF  CRUCIFORM  WING-BODY  ARRANGEMENT 


FIGURE  A-5B.  LOCATION  OF  WING-BODY  VORTICES  IN  PLANE  OF  TAIL 
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solutions  for  the  vortices  can  be  added  together  linearly,  the  analysis  for  the  <&  =  0  deg  plane  can 
be  utilized  for  the  $  =  45  deg  plane.  Referring  to  Figure  A-5A,  the  vortex  strengths  for  the  wing 
vortices  are  defined  by 


=  2  sin  «  cos  O 


(  s'  -  rM 


(A-25) 


=  -  Fj  =  2  V„  sin  «  cos  0 


'w 


/ 


Of  course  these  strengths  are  assumed  constant  from  the  wing  (Figure  A-5A)  to  the  tail  region 
shown  in  Figure  A-5B.  However,  both  f  or  y  and  h  or  z  change  for  each  of  the  vortices  in  going 
from  the  wings  to  the  tail  plane.  For  the  O  =  0  deg  case,  only  the  hs  changed. 

Referring  to  Figure  A-5B,  the  numbers  in  circles  refer  to  the  vortex  number  from  1  to  8, 
whereas  the  numbers  uncircled  refer  to  the  tail  fin  number.  In  analogy  to  the  <&  =  0  deg 
methodology,  the  impact  of  the  eight  vortices  on  tail  fins  1  and  4  is  required.  These  results  can  then 
be  doubled  to  account  for  fins  2  and  3.  Using  subscript  notation,  where  the  subscript  denotes  the 
vortex  number,  the  coordinates  of  the  wing-shed  vortices  with  respect  to  the  y-z  axis  system  are 


Zj  =  -  sin  0  +  tan  « 
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'  Tl'' 
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sin  (6^)  cos  0  {(c^)^/2 
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Yi  =  cos  $ 


(A-27) 
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(A-28) 

(A-29) 

Zj  =  Zj  +  2  cos  0  (A-30) 

Ya—yi  (A-31) 

Z4  =  Z3  (A-32) 


y4  =  Yi 


(A-33) 


Likewise,  the  coordinates  of  the  body  image  vortices  are  defined  by  (where  Equations  (A-13)  and 
(A- 14)  have  been  utilized) 


(A-34) 
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Yt  =  -  Yi 
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(A-39) 


(A-40) 


Yg  =  -  Y7  (A-41) 

Referring  back  to  Equations  (A- 16)  through  (A-23)  for  the  $  =  0  deg  roll  position,  note  that  the  lift 
on  the  right  tail  surface  was  computed  with  the  addition  of  four  vortices.  At  $  =  45  deg  roll,  there 
are  eight  vortices  that  act  on  the  two  tail  fins.  As  a  result,  the  wing  tail  interference  factor  for  tail 
one  is  (on  following  page) 
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Note  the  circulation  of  the  vortex  is  positive  in  a  counterclockwise  direction.  This  is  what 
determines  the  plus  or  minus  signs  in  Equation  (A-42).  Likewise  the  interference  factor  for  tail  4 
is  (on  following  page) 
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From  Equation  (10),  the  total  interference  lift  on  the  tail  surfaces  is 


'N. 


TtV) 


(A-44) 


It  is  interesting  to  note  that  for  zero  control  deflections,  the  term  given  by  Equation  (A-43) 

^T(V) 

should  be  independent  of  roll  position  This  means  linear  theory  gives  aerodynamics  of  a  wing- 
body-tail  configuration  that  are  independent  of  roll  position. 
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SUMMARY  OF  KEY  FINDINGS  FOR  ROLL  DEPENDENCE  WITH  RESPECT 
TO  LT  AND  SET 

The  whole  purpose  of  a  fairly  thorough  review  of  the  LT  and  SET  was  to  help  glean  the 
physics  involved  in  roll  orientation  at  low  angle  of  attack  in  order  to  help  structure  the  methodology 
that  will  incorporate  nonhnear  effects  at  larger  angles  of  attack.  As  such,  it  is  worthwhile  to  briefly 
summarize  these  key  results. 

a.  For  cruciform  wings  alone  or  a  wing-body  combination,  the  total  normal  force  is 
independent  of  roll. 

b .  For  a  planar  wing-body  combination  at  roll,  the  loading  on  the  windward  plane  panel 
is  greater  hy  an  equal  amount  to  that  on  the  leeward  plane  panel.  This  means  that 
if  one  were  trying  to  design  a  code  for  lateral  aerodynamics,  roll  dependence  of  each 
fin  planform  must  be  considered.  On  the  other  hand,  if  longitudinal  aerodynamics 
are  of  primary  interest,  the  total  normal  force  on  the  entire  wing  planform  can  be 
considered 

c.  For  a  cruciform  wing-body-tail  configuration  at  roll,  eight  vortices  are  shed  by  the 
wing-body  region,  which  adversely  affects  the  tail  lift.  This  is  as  opposed  to  four 
vortices  at  $  =  0  deg 

d.  The  planar  theory  developed  for  wing-tail  interference  can  be  used  to  approximate 
the  loss  of  lift  on  the  tails  at  $  =  45  deg 

e.  The  aerodynamics  of  a  cruciform  wing-body-tail  combination  with  zero  control 
deflections  are  independent  of  roll  position. 
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